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A new  methodology  is  developed  for  the  description  of  electronic  rearrange- 
ment in  atomic  and  molecular  collisions.  Using  the  eikonal  representation  of  the 
total  wavefunction,  time-dependent  equations  are  derived  for  the  electronic  densi- 
ties within  the  time-dependent  Hartree-Fock  approximation.  An  averaged  effective 
potential  which  ensures  time  reversal  invariance  is  used  to  describe  the  effect  of 
the  fast  electronic  transitions  on  the  slower  nuclear  motions.  Electron  transla- 
tion factors  (ETF)  are  introduced  to  eliminate  spurious  asymptotic  couplings,  and 
a local  ETF  is  incorporated  into  a basis  of  traveling  atomic  orbitals.  A refer- 
ence density  is  used  to  describe  local  electronic  relaxation  and  to  account  for  the 
time  propagation  of  fast  and  slow  motions,  and  is  shown  to  lead  to  an  efficient 
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integration  scheme.  Expressions  for  time-dependent  electronic  populations  and 
polarization  parameters  are  given.  Electronic  integrals  over  Gaussians  including 
ETFs  are  derived  to  extend  electronic  state  calculations  to  dynamical  phenomena. 

Results  of  the  method  are  in  good  agreement  with  experimental  data  for  charge 
transfer  integral  cross  sections  over  a projectile  energy  range  of  three  orders  of 
magnitude  in  the  proton-Hydrogen  atom  system.  The  more  demanding  calcula- 
tions of  integral  alignment,  state-to-state  integral  cross  sections,  and  differential 
cross  sections  are  found  to  agree  well  with  experimental  data  provided  care  is 
taken  to  include  ETFs  in  the  calculation  of  electronic  integrals  and  to  choose  the 
appropriate  effective  potential.  The  method  is  found  to  be  in  good  agreement 
with  experimental  data  for  the  calculation  of  charge  transfer  integral  cross  sec- 
tions and  state-to-state  integral  cross  sections  in  the  one-electron  heteronuclear 
Helium(2-t-)-Hydrogen  atom  system  and  in  the  two-electron  system.  Hydrogen 
atom -Hydrogen  atom. 

Time-dependent  electronic  populations  are  seen  to  oscillate  rapidly  in  the 
midst  of  collision  event.  In  particular,  multiple  exchanges  of  the  electron  are 
seen  to  occur  in  the  proton-Hydrogen  atom  system  at  low  collision  energies. 
The  concepts  and  results  derived  from  the  approach  provide  new  insight  into  the 
dynamics  of  nuclear  screening  and  electronic  rearrangement  in  atomic  collisions. 
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CHAPTER  1 
INTRODUCTION 


Electronic  rearrangement  in  ionic,  atomic,  and  molecular  collisions  is  a subject 
of  wide  interdisciplinary  interest.  It  deals  with  the  basic  concepts  and  methodol- 
ogy required  in  studies  ranging  from  detailed  atom-atom  collision  phenomena  to 
phenomena  in  very  complex  systems.  For  example,  electron  transfer  in  photoex- 
cited  proteins  and  protein-enzyme  interaction  is  of  current  interest  in  biological 
and  biochemical  research.  Electron  transfer  in  atom-surface  and  molecule-surface 
interaction  is  an  area  of  active  research  in  material  science  and  electron  transfer 
in  electrochemical  cells  is  of  interest  in  current  chemical  engineering  research. 
Thus,  it  is  of  far  reaching  importance  to  gain  a fundamental  understanding  of 
the  most  basic  time -dependent  aspects  of  electron  transfer  phenomena  in  simple 
systems  where  detailed  studies  can  be  carried  out.  This  thesis  investigates  these 
time-dependent  phenomena  in  terms  of  a formalism  developed  to  consider  many 
active  electrons. 


Overview  of  Theoretical  Methods 

While  the  formalism  developed  is  applicable  to  system  involving  many  atomic 
cores  and  electrons,  in  this  thesis  we  consider  the  most  basic  of  collision  in  which 
there  are  only  two  cores  and  one  or  two  active  electrons.  The  earliest  theoretical 


1 


2 


treatment  considered  in  this  thesis  is  the  Perturbed  Stationary  States  treatment 
which  dealt  specifically  with  two  identical  cores  and  only  one  electron  [1,  2]. 
Others  have  returned  to  this  formalism  more  recently  [3,  4].  In  this  treatment  the 
two  cores  are  taken  to  be  essentially  at  rest  with  the  motion  of  the  cores  taken 
to  be  a perturbation  which  induces  electronic  transitions.  For  example,  the  total 
wave  function  of  the  system  with  identical  nuclei  is  written 

= X+{p,R)F+{R)  + X-{p,R)F-{R)  (1.1) 

where  p is  the  electron  coordinate,  R is  the  intercore  separation  and  X'^{p,R) 
and  X~{p,R)  are  symmetric  and  antisymmetric  in  the  nuclear  coordinates  with 
coefficients  F~^{R)  and  F~{R),  respectively.  The  total  wave  function  is  substi- 
tuted into  the  time-dependent  Schrddinger  equation  and  the  solutions  are  required 
to  have  this  asymptotic  form 

F^{R)  ~ 4=[exp(titn  • R)  -t-  R~^f^{e)  exp(ikR)]  (1.2) 

V2 

where  k is  the  magnitude  of  the  initial  wave  vector  of  the  relative  motion,  n 
is  a unit  vector  in  the  direction  of  the  final  relative  motion,  0 is  the  scattering 
angle  and  /^(0)  are  the  scattering  amplitudes.  The  integral  charge  transfer  cross 
section  is  then  expressed  as 

^ = - r?s\nedQd^ 

0 0 


(1.3) 
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While  this  theory  obtains  good  results  for  integral  cross  sections  at  low  energies, 
its  assumption  breaks  down  at  higher  energies  and  the  inclusion  of  more  than  two 
total  states  for  the  wave  function  complicates  the  theory.  Perturbed  stationary 
states  also  fail  to  give  time  dependent  information  about  the  collision. 

Another  approach  to  the  dynamics  of  collisional  systems  has  been  the  direct 
numerical  integration  of  the  Schrodinger  equation  within  the  impact  parameter 
approximation  [5].  In  this  reference,  the  collision  system  has  been  assumed  to 
be  homonuclear  and  occurs  in  the  x-z  plane.  The  time-dependent  Schrodinger 
equation  is  then  written  in  cylindrical  coordinates  as 

1 d9 

idrp{p,z,<p,t)ldt  = + V{p,z,t)  - •^Ly]ip{p,z,ip,t)  (1.4) 

where  ‘ip{p,z,(p,t)  is  the  electronic  wave  function,  Ly  is  the  y-component  of  the 
angular  momentum  operator,  6 is  the  angle  of  rotation  between  the  space-fixed 
and  body-fixed  axes  in  the  x-z  plane,  and  the  Coulomb  potential  V is  given  by 

V(p,z,t)  = + |z  + - {p'‘  + (z  - 5H(<)1-)-’  (1-5) 

where  R{t)  is  the  intemuclear  separation.  Expansion  of  the  wave  function  in  the 
eigenfunctions  of  the  z component  of  the  angular  momentum  reduces  the  problem 
to  one  of  two  spatial  variables.  After  the  initial  wave  function  is  specified  on  a 
grid  in  the  p-z  plane,  the  numerical  integration  is  carried  out.  In  this  reference,  the 
nuclear  trajectory  is  assumed  to  be  a straight  line  with  constant  velocity.  While 
this  method  can  give  electronic  densities  as  a function  of  time,  it  does  not  give 
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insight  into  other  physical  properties.  The  restriction  of  straight  line  trajectories 
hinders  the  determination  of  differential  cross  sections  and  integral  cross  sections 
at  low  collision  energies.  At  higher  collision  energies,  the  method  gives  good 
total  integral  cross  sections  and  state-to-state  integral  cross  sections.  Finally,  the 
method  is  very  computationally  intensive  so  that  generalizing  it  to  larger  systems 
may  not  be  practical. 

By  far  the  most  popular  method  for  calculating  the  dynamic  of  collisional 
system  is  that  of  close-coupling.  Two  recent  reviews  give  a good  overview  of 
this  work  [6,  7].  In  this  method,  the  time  dependent  electronic  wave  function 
is  expanded  in  a basis  of  atomic  states,  molecular  states,  pseudo-states,  or  a 
combination  of  such  states 

N 

= '^ak{t)4^k{r,t)  (1.6) 

k=i 

where  the  expansion  coefficients  ak{t)  are  complex  and  time  dependent.  The 
expansion  coefficients  are  propagated  through  time  using  the  time-dependent 
Schrodinger  equation  along  a chosen  nuclear  trajectory.  In  these  calculations  the 
trajectory  choice  is  generally  the  straight  line  trajectory  with  constant  velocity. 
After  the  completion  of  the  nuclear  trajectory  the  transition  probability  to  a final 
state  / is  given  by 


P,^f{v,b)  = |a/(-|-oo)|2 


(1.7) 
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where  b is  the  impact  parameter  and  the  state-to-state  integral  cross  section  to 
this  state  is 


The  choice  of  basis  functions  is  what  distinguishes  the  various  close-coupled 
calculations. 

The  earliest  choice  of  basis  functions  was  the  hydrogenic  functions  [8-10]. 
All  of  these  studies  included  electron  translation  factors  and  hydrogenic  orbitals 
on  each  center.  The  role  of  electron  translation  factors  will  be  discussed  in  chapter 
3.  Various  types  of  atomic  orbitals  have  been  used  in  calculations  involving  two 
centers.  Some  of  these  consider  atomic  orbitals  appropriate  to  the  individual 
centers  only  and  are  called  AO  approaches  [11-13].  Others  include  not  only 
atomic  orbitals  appropriate  to  the  individual  centers  but  also  those  appropriate 
to  the  charge  of  the  united  atom.  These  are  called  AO-i-  approaches  [14-16]. 
Another  common  choice  of  basis  functions  is  a set  of  molecular  orbitals  [17-20]. 
It  is  more  difficult  to  decide  the  proper  choice  of  the  electron  translation  factor  in 
this  basis.  It  is  also  necessary  to  match  the  molecular  orbital  basis  to  an  atomic 
orbital  basis  asymptotically  to  obtain  cross  section  information.  Pseudo-states 
have  been  included  in  the  expansion  of  atomic  orbitals  including  Sturmians  and 
Hylleraas  functions  [21-23].  Finally,  there  has  been  an  attempt  at  improving  upon 
the  trajectory  choice  through  the  so-called  common  trajectory  method  [24,  25]. 
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While  close-coupled  calculations  give  good  agreement  with  experimental 
results  at  high  energies,  the  generalization  of  this  method  to  more  than  two  centers 
is  quite  challenging.  Also  the  calculation  of  time  dependent  quantities  of  physical 
interest  has  not  been  done  within  the  close-coupled  method. 

The  present  method  combines  and  builds  upon  methods  developed  over  a 
number  of  years  using  eikonal  procedures  and  time-dependent  Hartree-Fock  states 
[26-30].  The  eikonal  representation  of  the  total  wave  function  of  a system  and 
its  relationship  to  the  eikonal  approximation  have  been  developed  particularly  in 
the  context  of  fast  and  slow  degrees  of  freedom.  Time-dependent  Hartree-Fock 
has  been  the  subject  of  much  investigation  particularly  into  the  circumstances 
in  which  time-dependent  Hartree-Fock  must  be  improved.  A general  variational 
method  based  on  a transition  functional  has  provided  an  improvement  scheme 
and  has  been  applied  to  the  Helium(2-t-)-Helium  collisional  system  [31].  The 
development  of  the  present  method  was  initially  concerned  with  the  equations  for 
the  time  propagation  of  the  electronic  densities  within  time -dependent  Hartree- 
Fock  theory  and  with  the  coupling  of  electronic  and  nuclear  motions  [32-34]. 
In  this  method,  we  use  the  equations  for  the  time  propagation  of  the  electronic 
density  operator 


(1.9) 


where  p is  the  electronic  density  operator  and  F is  the  Fock  operator.  The 
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choice  of  propagating  the  electronic  densities  is  motivated  by  the  fact  that  the 
electronic  density  is  more  fundamental  to  the  physics  of  collisional  systems  than 
the  expansion  of  the  wave  function  in  molecular  orbitals.  The  electronic  density 
also  has  the  advantage  of  being  closely  related  to  physical  quantities.  The  method 
has  been  further  developed  to  consider  such  physical  quantities  as  the  time- 
evolution  of  electronic  polarization  and  electronic  populations  and  to  calculate 
collisional  polarization  parameters  [35,  36]. 

An  approach  of  interest  related  to  this  thesis  is  based  on  the  time-dependent 
variational  principle  and  coherent  states  [37,  38].  It  has  recently  been  applied  to 
the  electron-nuclear  dynamics  of  atom-atom  collisions  [39].  This  theory  is  allied 
to  the  present  method  in  that  it  is  also  a time-dependent  Hartree-Fock  method  in 
which  the  nuclear  and  electronic  degrees  of  freedom  are  coupled.  The  formalism 
is  developed  in  a coherent  states  representation  so  that  wave  function  takes  the 
form 


creation  and  annihilation  operators,  respectively,  and  the  complex  parameters 
Zph  called  the  Thouless  parameters  define  the  state  of  the  system.  The  time- 
dependent  Hartree-Fock  equations  are  derived  as  a special  case  of  a more  general 
time-dependent  variational  principle  and  the  Thouless  parameters  are  propagated 


(1.10) 


where  |0)  is  a single-determinantal  reference  wave  function,  6p  and  bh  are  the 
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through  time.  The  coupling  of  the  nuclear  and  electronic  motion  is  accomplished 
through  the  use  of  a time  reversal  invariant  generalization  of  the  Ehrenfest 
potential  to  account  for  the  effect  of  the  electronic  motion  on  the  nuclear  motion. 

Outline  of  the  Thesis 

Chapter  2 deals  with  the  coupling  of  electronic  and  nuclear  degrees  of  freedom 
within  the  eikonal  representation.  The  electronic  degrees  of  freedom  are  treated 
quantally  and  the  potential  the  cores  experience  due  to  the  electrons  is  formulated 
in  a time  reversal  invariant  form.  Definitions  of  the  types  of  electronic  potentials 
used  throughout  the  thesis  are  given.  Chapter  3 presents  the  derivation  of 
the  time-dependent  Hartree-Fock  equations  for  the  electronic  density  matrices. 
The  derivation  includes  a detailed  exposition  of  the  formal  role  of  the  electron 
translation  factor  in  the  context  of  the  present  approach.  Chapter  4 develops  a 
method  for  the  solution  of  these  equations  by  a temporal  linearization  technique. 
A solution  of  the  linearized  equations  by  an  exponential  expansion  is  presented, 
and  practical  criteria  for  the  efficient  integration  of  the  equations  are  advanced. 
Chapter  5 discusses  time-dependent  properties  for  the  more  detailed  description 
of  collisional  dynamics.  Particular  attention  is  given  to  measures  of  electronic 
population  and  polarization  parameters. 

Chapter  6 sets  out  the  procedure  for  the  integration  of  the  linearized  time- 
dependent  Hartree-Fock  equations  along  with  the  integration  of  the  nuclear 
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equations  of  motion.  A scheme  is  developed  which  allows  for  the  rapid  inte- 
gration of  these  equations  without  a loss  of  accuracy.  Chapter  7 deals  with  our 
basis  set  choices  and  their  effects  on  our  calculations.  The  role  of  electron  trans- 
lation factors  at  differing  core  velocities  is  considered.  A method  for  the  inclusion 
of  electron  translation  factors  in  overlap  and  Hamiltonian  integrals  involving  two 
centers  is  developed  for  Gaussian-type  orbitals. 

Chapter  8 discusses  the  results  of  the  application  of  our  method  to  the  proton- 
Hydrogen  collisional  systems.  Results  are  presented  over  a range  of  projectile 
energies  from  a few  eV  to  a few  keV  for  charge  transfer  integral  cross  section, 
state-to-state  integral  cross  sections,  differential  cross  sections,  and  integral  align- 
ment. Time  evolution  of  electronic  populations  and  alignment  parameters  is  also 
presented.  Chapter  9 contains  results  for  charge  transfer  integral  cross  section 
and  state-to-state  integral  cross  sections  for  the  Helium(2+)-Hydrogen  collisional 
system.  Again,  time  evolution  of  electronic  populations  and  alignment  parame- 
ters are  shown.  Chapter  10  considers  the  two  active  electron  Hydrogen-Hydrogen 
collisional  system.  Integral  state-to-state  cross  sections  and  time  evolution  of 
electronic  populations  are  presented.  Chapter  1 1 discusses  the  method  and  de- 
scribes possible  improvement  to  the  methods  and  future  directions.  Appendices 
A and  B present  the  details  of  the  calculation  of  overlap  and  Hamiltonian  matrix 
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elements  in  the  Hydrogenic  and  travelling  Gaussian  bases,  respectively.  Appendix 
C presents  a flow  diagram  of  the  time-dependent  density  matrix  code. 


CHAPTER  2 

COUPLING  OF  ELECTRONIC  AND  NUCLEAR  MOTIONS 


The  Hamiltonian  operator  and  the  total  wavefunction  of  a collisional  system 
can  be  written  in  atomic  units  as 

^ = + (2.,) 

^'(X,R)  = x(X,R)exp[i5(R)] 

where  X and  R are  the  active  electronic  and  core  coordinates.  The  second  term 
in  the  Hamiltonian  describes  the  active  electronic  system  for  fixed  cores.  The 
function  5 is  real,  while  x is  complex-valued;  in  our  treatment  we  choose  S 
for  convenience  to  be  an  action  integral  approximating  the  quantal  phase  of  the 
wavefunction;  the  phase  of  x can  compensate  in  principle  for  deviations  from  the 
correct  phase.  We  refer  to  this  decomposition  of  the  wavefunction  as  the  eikonal 
representation. 

For  total  energy  E,  we  place  the  wavefunction  in  the  time-independent 
Schrodinger  equation  to  obtain 


(2.2) 


which  is  valid  for  any  S.  A specific  equation  for  S follows  by  projecting  equation 
2.2  on  X.  and  taking  the  real  part  of  the  result,  which  gives  [26,  29] 


2M^dR^  ^ 

= (x(R)|^r|x(R))/{x(R)Ix(R))  + 


(2.3) 
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where  and  contain  5x/5R  and  respectively,  and  the  first 

term  can  be  identified  as  the  Ehrenfest  potential  for  the  core  motions. 

Defining  the  momentum 

Pr  = dSIdK  (2.4) 


trajectories  of  the  cores  follow,  introducing  the  time  t,  from 
H,u  = Pr  • Pr/(2M)  + V,„(Pr,  R) 

dRidt^  a//,„/aPR 


dPKldt=  -dHguldn 
for  initial  (PR^,Ra)  and  give 


R 

5(R)  = 5(Ra)  + J dR  PK 

Ra 


(2.5) 


(2.6) 


We  assume  in  this  treatment  that  the  electronic  states  can  be  obtained  from 
the  time-dependent  Hartree-Fock  (TDHF)  approximation.  This  includes  electron 
correlation  (beyond  the  time-independent  HF  approximation),  and  electronically 
diabatic  effects,  insofar  as  it  allows  for  electronic  transitions  induced  by  the  core 
motions.  We  further  assume  that  we  can  obtain  the  effective  potential  within  an 
eikonal  (short  wavelength)  approximation.  We,  therefore,  neglect  and 
in  equation  2.3  and  write  the  effective  potential 


Vgu=  tr[p 


(Me) 


2 


(2.7) 
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however,  as  the  Fock  operator,  F,  the  core  Hamiltonian  He  and  the  density 
operator  p are  each  functions  of  time,  this  Ehrenfest  potential  is  neither  manifestly 
nor  generally  time  reversal  invariant  for  state-to-state  transitions.  An  averaged 
effective  electronic  potential  can  be  written  in  a form  invariant  under  time 
reversal 

Vqu  = tr  (^  + (2.8) 

li 

where  p denotes  the  initial  electronic  state  of  the  system  and  the  are  appro- 
priately chosen  weights.  This  potential  will  be  time  reversal  invariant  if  all  the 
appropriate  active  channels  are  considered.  Criteria  for  the  choice  of  appropriate 
weight  functions  include 

f 

1.  All  open  channels  in  the  case  where  only  open  channels  participate  in  a given 
collision. 

2.  All  resonant  virtual  excitation  channels  which  predominate  in  a given  colli- 
sion. 

In  this  work  we  have  employed  a number  of  approximate  potentials  for  the 
motions  of  the  cores.  The  first  approximation,  known  as  the  Impact  Parameter 
Approximation  , consists  of  allowing  the  core-core  potential  and  the  core-electron 
potential  to  be  zero  (ZP).  The  result  is  a rectilinear  trajectory  for  the  cores  with 
constant  velocity.  This  approximation  is  best  used  at  high  collision  energies  and 
is  used  for  the  purpose  of  comparison  to  other  results  in  the  literature. 
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A second  approximate  potential  we  will  employ  we  call  the  Coulomb  Potential 
(CP).  It  consists  of  keeping  the  core  repulsion  while  setting  the  core-electron 
potential  to  zero. 

V—  Vcore  — CaCbIRab  (2-9) 

Here  (at  is  the  effective  charge  on  center  N.  This  approximation  will  be  used  to 
compare  with  the  Impact  Parameter  Approximation  and  other  more  exact  forms 
for  the  core-electron  potential.  One  would  expect  that  the  Coulomb  potential 
would  be  most  successful  at  fairly  high  collision  energies. 

A third  approximate  potential  we  will  call  the  Screened  Coulomb  Potential 
(SCP).  This  choice  of  potential  consists  of  keeping  the  core  repulsion  while 
choosing  the  weights  from  equation  2.8  such  that  the  core  screening  by  the 
valence  electrons  for  the  ground  states  in  a given  system  is  included  while  all 
other  electronic  states  are  ignored.  This  potential  is  especially  successful  when 
lower  collision  energies  are  involved  and  particularly  for  low  collision  energy, 
resonant  transfer  phenomena.  It  will  be  shown  to  give  results  in  good  agreement 
with  sensitive  experimental  quantities  such  as  differential  cross  sections. 

A fourth  approximate  potential  will  be  the  Averaged  Effective  Potential  (AEP). 
This  choice  of  potential  involves  retaining  the  core  repulsion  while  choosing 
weights  Wfi  from  equation  2.8  such  that  all  open  channels  for  a given  collision  are 
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equally  weighted.  This  potential  should  be  successful  at  intermediate  and  lower 
collision  energies  and  particularly  for  nonresonant  transfer  collisional  phenomena. 


CHAPTER  3 

THE  TIME-DEPENDENT  HARTREE-FOCK  EQUATIONS 


The  inclusion  of  time  dependence  in  the  theory  of  ionic,  atomic,  and  molecular 
collisions  entails  a number  of  challenges.  Unlike  the  time-independent  theory, 
time-dependent  Hartree-Fock  needs  to  consider  a large  range  of  core  geometries 
and  nonequilibrium  electronic  configurations.  Further  the  integration  of  the  time- 
dependent  Hartree-Fock  equations  must  allow  for  rapid  variations  in  electronic 
densities  around  a given  core.  It  would  also  be  helpful  if  these  equations  offered 
a unique  physical  insight  into  the  dynamics  of  collisional  systems  involving  both 
fast  and  slow  degrees  of  freedom. 

A collisional  system  including  ions,  atoms,  and  molecules  will  be  character- 
ized by  a mathematical  description  of  the  core  and  electronic  degrees  of  freedom. 
In  the  present  method  the  core  degrees  of  freedom  will  be  treated  classically  (see 
Chapter  2).  The  electronic  degrees  of  freedom  will  be  described  quantally  through 
the  time-dependent  Hartree-Fock  equation.  One  further  needs  to  decide  whether 
one  wishes  to  treat  the  electronic  degrees  of  freedom  by  the  propagation  of  their 
wave  function  through  time  or  an  alternative  quantity.  In  the  present  method 
we  opt  to  propagate  the  electronic  density  matrices  through  time  as  opposed  to 
propagating  the  electronic  wave  function  or  some  other  amplitude. 
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The  time-dependent  Hartree-Fock  equations  for  the  electronic  densities  can 
be  written  in  operator  notation  with  ft  = 1, 

i'p  = FV  - (3.1) 

where  the  dot  denotes  differentiation  with  respect  to  time, 

= (3-2) 

occ  i 

is  the  electronic  density  operator  for  electron  spin  7 and  <f>-  are  the  molecular 
orbitals  which  compose  the  time-dependent  Hartree-Fock  state,  and  the  Fock 
operator  is  defined  by 

F-^  = H.  + G^p^p-^']  (3.3) 

with  the  single  electron  Hamiltonian  H^,  and  the  electron-electron  interaction  G"*^. 

We  expand  the  molecular  orbitals  as  linear  combinations  of  traveling  atomic 
orbitals  (TAOs).  We  implement  this  expansion  by  writing  the  molecular  orbitals 

(3.4) 

Here  we  have  introduced  the  local  version  of  the  electron  translation  factor  (ETF) 
[40,  41] 

t 

Tm{rj,t)  - exp[ime{vm  J 

to 


(3.5) 
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where  the  local  velocity  vector  for  core  m is  the  position  of  the  jth  electron 
with  respect  to  a space-fixed  origin  is  fj,  the  position  of  the  jth  electron  with 
respect  to  core  m is  fm],  the  mass  of  the  electron  is  me,  an  atomic  orbital 
centered  at  m is  X/x,  and  the  expansion  coefficients  are  c^-.  The  role  of  the  electron 
translation  factor  is  to  cancel  spurious  couplings  that  arise  from  the  kinetic  energy 
operator  acting  on  atomic  states  whose  origins  are  functions  of  time.  Here  we 
will  consider  only  situations  in  which  the  asymptotic  accelerations  of  the  cores 
are  zero.  This  construction  can  be  modified  to  include  asymptotically  accelerating 
particles.  The  single-electron  Hamiltonian  can  be  written  as  a sum  of  kinetic  and 
potential  energy  operators 

H = f+V  (3.6) 

where  the  kinetic  energy  operator  is  defined  by 

(3.7) 

Consider  the  action  of  the  time  derivative  on  an  atomic  state  with  an  origin  that 
is  a function  of  time  in  the  asymptotic  region. 

X/x(^mj)  = *[(^){^^}  + X/iO>ni)  (3.8) 

where  the  brackets  indicate  which  quantities  are  held  constant  while  the  partial 
time  derivative  is  taken.  Note  the  presence  of  the  gradient  with  respect  to  the  core 
coordinate  Rm  on  the  right  hand  side  of  equation  3.3.  This  is  often  referred  to 
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as  the  spurious  gradient  coupling  which  causes  the  mixing  of  pure  atomic  states 
on  the  core  m.  Introducing  the  TAO 


(3.9) 


applying  the  time  derivative  operator  to  a TAO  we  obtain 

'1 
^dt 


Ur„t)  = ^ ■ VjiJ 


dt 


(3.10) 


where  we  employ  the  fact  that  fj  is  independent  of  time.  Continuing  in  more 
detail 

*(^){f_,}l^mX>)  = ^ TmlXfi)  + 


p IX/i) 

,2  (3-H) 

T\Tmx,)  = - idm  • V,-  + r^r]|x,)  = 

[-^Y-Trn  + iTmVm  ' + r^rllX/.) 

Here  one  can  readily  see  that  in  the  limit  where  the  acceleration  can  be  neglected 
the  first  and  third  terms  from  the  time  derivative  cancel  the  first  and  second  term 
in  the  kinetic  energy,  respectively.  This  leaves  the  usual  time-dependent  Hartree- 
Fock  equation  in  the  ^-basis,  but  with  an  altered  kinetic  energy  operator 


irk,)  = r„7T-'|f„) 


(3.12) 


Hence,  the  modified  Fock-like  operator  can  be  defined  by 


(3.13) 
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which  is  not  generally  a Hermitian  operator.  We  may  now  write  the  time- 
dependent  Hartree-Fock  equation  in  the  TAO  basis 

f:}p^  - = isp^s^ 


(3.14) 


s"^F].p^  - p^f:[^s-^  = ?p^ 

where  bold  type  indicates  matrices  in  the  TAO  basis,  S“^  is  the  inverse  of  the 
overlap  matrix,  Fj  is  the  modified  Fock-like  matrix  and  the  density  matrix  is 
given  by  its  usual  definition  in  terms  of  the  expansion  coefficients 


c«)=  E (3.15) 

occ  i 

In  the  next  chapter  we  will  present  a method  of  solution  of  these  equations  for 
the  propagation  of  the  electronic  density  by  the  method  of  temporal  linearization. 


CHAPTER  4 

TEMPORAL  LINEARIZATION  OF  THE 
TIME-DEPENDENT  HARTREE-FOCK  EQUATIONS 

In  the  previous  chapter,  we  derived  the  Hartree-Fock  equations  for  the  elec- 
tronic densities  in  matrix  form  for  a basis  of  traveling  atomic  orbitals.  In  this 
chapter  we  will  present  a method  for  solution  of  these  equations  employing  the 
technique  of  temporal  linearization.  The  linearization  procedure  can  be  presented 
in  operator  notation.  Let  us  begin  by  writing  the  electronic  density  operator  in 
the  time  interval  between  to  and  t as 

P’(0  = P<">"(0 + ■?’(<)  (4.1) 

Further  define  the  reference  density  to  satisfy  the  equation 

(4.2) 

Here  we  need  to  make  some  remarks  concerning  the  modified  Fock-like  operator 
which  can  be  written 

F}[p\ p)  = TV  + K -f  G[p\ p)  (4.3) 

where  TV  is  the  modified  kinetic  energy  term,  Vc  is  the  core  potential,  and 
P ] is  the  electron-electron  interactions  among  active  electrons.  The  choice 
of  core  potential  and  definition  of  active  electrons  will  vary  from  system  to  system 
and  will  also  depend  on  the  collision  energies  to  be  considered. 
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Given  these  definitions,  we  can  then  derive  an  equation  for  the  change  in 
the  density  matrix  a‘’{i)  and  keep  terms  only  to  first  order  in  this  change  in  the 
density  matrix. 

ia\t) 

(4.4) 

where 


AH^t)  = frit)  + Vc{t)  - [frito)  + K(<o)] 


(4.5) 

AG[p(0)>,p(0)V]  = G[p(0)7(<),p(0)7'(^)]  _ G[^(0)T(<o),^(0)7'(^^)| 

We  can  now  see  that  the  reference  density  that  we  have  defined  has  the  effect 
of  separating  the  electronic  density’s  evolution  arising  from  the  TDHF  equation 
with  the  cores  fixed,  from  the  electronic  density’s  evolution  that  arises  from  the 
effect  of  the  movement  of  the  cores.  More  specifically  the  reference  density 
evolves  as  though  the  cores  are  fixed  while  the  change  in  the  density  a'*{t)  corrects 
for  the  effect  of  the  movement  of  the  cores. 

These  equations  move  easily  into  the  TAG  basis  matrix  form,  with  equation 
4.2  for  the  reference  density  becoming 


,p(0)7(f)  = w^(^o)P(°)^(0  - (4.6) 


where 

= S-^F^ 


(4.7) 
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Continuing  along  this  path,  equations  4.4  and  4.5  for  the  change  of  the  density 
may  be  written  in  matrix  form 


= W^(to)Q^(/)  - + D^(0 

D^(t)=  AWWT'(t)pWT'(t)  - 
where  in  the  driving  term  D'*'(t)  we  define 


(4.8) 


W(°h  = + G[P(°)\  p(o)Vj j (4  9) 


The  matrix  equations  for  the  electronic  densities  are  a set  of  linear  first- 
order  differential  equations.  The  solution  to  these  equations  is  problematic  as  the 
change  in  the  density,  as  well  as  the  reference  density,  may  oscillate  significantly. 
Therefore,  we  choose  to  solve  these  equations  using  an  exponential  expansion.  In 
this  procedure,  the  change  in  the  density  matrix  is  expressed  as 

Q'^(<)=  exp(-iW'^t)Q'(t)exp(iW^^t)  (4.10) 

Substituting  equation  4.10  into  equation  4.8,  we  obtain 

eexp(-zW^t)Q'(t)exp(zW^tt)=  D^(t)  (4.11) 

Returning  to  Q'’'(f),  its  expression  at  time  ti  is 
Q^(^a)=  -iexp(-iW'^t,)x 

[ j exp{iW'ft')D\t')exp{-iW‘^U')di']  exp(tW^^t,) 

<0 


(4.12) 
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Introducing  the  biorthogonal  transformations  which  diagonalize  and 
[42] 

= Lw^L^ 

(4.13) 


= Lw^L^ 

We  may  now  write  the  matrix  elements  of  the  change  of  the  density 


jpak 


kv 


(4.14) 


where 


rjp<TJfc(^i’^o)  = J exp[-i{w]  - ivl*){h  - i')]D],„{t')di'  (4.15) 

<0 

where  the  w'J  are  complex  eigenvalues  of  the  matrix  W'>'.  Hence,  we  have 
replaced  the  single  integral  over  time,  from  the  initial  time  to  the  final  time,  by  a 
series  of  integrations  from  to  to  ti  where  to  is  the  initial  time  and  ti  is  the  final  time 
corresponding  to  the  chosen  ending  separation  between  the  cores.  The  advantage 
of  such  an  approach  is  that  the  problem  of  rapid  variations  in  the  density  matrix 
is  avoided  by  accounting  for  such  variations  explicitly  in  the  diagonalization  of 
the  W')'  and  matrices. 

In  the  case  of  two  or  more  active  electrons  a tolerance  for  iterative  solution 
of  the  TDHF  equations  will  be  considered.  If  the  ratio  of  the  magnitude  of  the 
change  of  the  density  matrix  to  the  magnitude  of  the  density  matrix  itself  is 
greater  than  a given  tolerance 

IIQ"(<i)ll 


^iter 


(4.16) 
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we  return  to  time  to  substituting 

K+i{k)=  PlitoHfQlih) 

o</<  1 


(4.17) 


where  the  subscript  indicates  the  number  of  iterations,  to  define  a new  starting 
point  for  the  reference  density.  The  choice  of  the  fractional  weight  / is  made 
to  optimize  the  convergence  of  the  iterative  procedure.  The  procedure  can  be 
portrayed  schematically 


p;e.(<o)=  p;„(<o)+/q:„(;i) 


else  accept 


(4.18) 


Here  the  first  line  is  repeated  until  the  iteration  criterion  is  satisfied. 

The  above  separation  of  the  density  matrix  into  a reference  density  and  a 
change  in  the  density  matrix  leads  naturally  to  a set  of  criteria  for  the  appropriate 
choice  of  to  and  ti.  A time  step 


At  = — to 


(4.19) 


will  be  deemed  acceptable  provided  the  ratio  of  the  magnitude  of  the  change  of 
the  density  matrix  to  the  magnitude  of  the  density  matrix  itself  is  between  given 
tolerances, 

Slower  < ||P7(^j)||  ^higher  (4.20) 

Here  the  lower  tolerance  tiowtr  allows  the  method  to  take  larger  time  steps  in 
any  region  in  which  the  evolution  of  the  density  matrix  is  essentially  equal  to 
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the  evolution  of  the  reference  density.  The  higher  tolerance  Chigher  restricts  the 
method  to  take  smaller  time  steps  in  any  region  in  which  the  change  in  the  density 
matrix  is  of  great  importance,  thereby  keeping  the  error  in  the  procedure  within 
a prescribed  level.  The  present  method  is  therefore  able  to  automatically  adjust 
its  time  step  size  in  order  to  maximize  the  performance  of  the  method  in  terms 
of  taking  the  minimum  number  of  possible  steps  for  the  maximum  accuracy  for 
those  steps.  In  practice  the  choice  of  tolerances  will  be  shown  to  be  rather  flexible 
in  achieving  the  desired  accuracy  in  an  acceptable  amount  of  computing  time. 

In  the  case  of  one  active  electron,  we  simply  replace  the  Fock  operator  with 
the  full  Hamiltonian  for  the  active  electron.  So  we  write,  for  a given  spin, 

W = S-^Ht  (4.21) 

the  remainder  of  the  development  proceeds  as  above.  There  is  no  need  for 
iterative  propagation  in  the  one  active  electron  limit.  This  method  will  be  used 
to  investigate  electron  transfer  and  excitation  in  systems  with  one  active  electron. 
We  will  use  this  method  to  investigate  polarization  and  alignment  effects  in 
collisions.  The  effect  of  electron  translation  factors  and  differing  basis  sets  will  be 
investigated.  We  will  also  explore  candidates  for  the  description  of  time-dependent 
phenomena  during  atom-atom  and  atom-ion  collisions  which  will  serve  as  the 
topic  of  the  next  chapter. 


CHAPTER  5 

TIME-DEPENDENT  ELECTRONIC 
ALIGNMENT  AND  POPULATION  ANALYSES 


It  is  of  interest  to  describe  in  some  way  what  is  occurring  during  collisional 
electronic  excitation  and  rearrangement.  Among  the  possible  quantities  of  interest 
are  the  alignment  and  orientation  parameters  of  the  collisional  system  and  the 
electronic  populations  of  the  atomic  and  molecular  orbitals  of  the  system.  While 
a number  of  descriptions  appropriate  to  time-independent  calculations  are  available 
in  the  literature,  we  will  propose  a number  of  time-dependent  extensions  of  some 
of  these  descriptions.  We  will  evaluate  these  extension  by  applying  them  to 
a number  of  differing  collisional  systems  for  a variety  of  collision  energies. 
Further,  we  will  discuss  the  applicability  of  these  measures  of  electronic  alignment 
and  population  with  particular  emphasis  on  physical  insights  gained  from  time- 
dependent  studies. 

In  this  study,  we  will  always  consider  the  plane  of  the  collision  to  be  the  x-z 
plane  with  the  z-direction  being  the  incoming  direction  of  the  projectile  and  the 
x-direction  being  the  direction  of  the  impact  parameter,  so  that  the  y-direction  is 
always  perpendicular  to  the  plane.  A few  other  investigators  choose  the  plane  of 
the  collision  to  be  the  x-y  plane  [43,  44],  but  the  consensus  of  the  field  seems  to 
agree  with  our  calculational  geometry  [45,  46]. 
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projectile 


28 


29 


Having  chosen  this  particular  geometry,  we  may  embark  upon  the  definition 
of  orientation  and  alignment  parameters.  In  general  one  may  describe  the  orbital 
polarization  of  a system  by  examining  the  anisotropy  of  the  orbital  angular 
momentum  of  the  electrons.  A state  of  given  angular  momentum  L is  said  to 
be  [47] 

1.  Isotropic,  if  all  its  magnetic  substates  are  equally  populated. 

2.  Oriented,  if  the  magnetic  substates  have  differing  populations. 

3.  Aligned,  if  substates  with  magnetic  quantum  numbers  of  equal  magnitude  but 
opposite  signs  are  equally  populated. 


To  illustrate  these  definitions  let  us  consider  a system  with  total  angular  momentum 
J.  We  may  introduce  the  spherical  irreducible  tensor  operators  where  we 
will  consider  n=l,2  and  -n  < q < n 


— j 

Jq  — Jz 
7(1)  _ J ^ 7 


7(2)  _ _1 


jr  = -^(3J2  _ j2) 


J^i  — ^J±{2Jz  ± 1) 


j{2)  _ j.  j2 

J±2  - -'J- 


where  we  define  in  Cartesian  coordinates 


— Jx  i 


(5.1) 


(5.2) 
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The  action  of  the  operators  is  defined  by 

|J.,  if  ] = 

= VHk  + 1)  - til  ± 1) 


(5.3) 


where  the  square  brackets  indicate  the  commutator.  Let  us  now  define  the 
orientation  vector  and  the  alignment  tensor  at  time  t, 


0<‘V,()  = 


1 

^/AJ+T) 


)t 


(5.4) 


(5.5) 

where  the  angle  brackets  indicate  the  expectation  value  with  respect  to  p{t)  for 
a given  J.  These  quantities  can  be  expressed  in  terms  of  the  density  matrices 
and  hence  may  be  calculated  as  functions  of  time  in  our  method.  Examining  the 
expression  for  the  expectation  value  of  the  spherical  irreducible  tensor  operators 
we  may  write 

( (5.6) 

where  the  subscript  m indicates  the  core  to  which  the  tensor  will  be  referred.  This 
expression  can  be  written  in  terms  of  the  expansion  in  TAOs 

fir 
fiv 


(5.1) 
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where  we  have  employed  the  cyclic  property  of  the  trace.  We  further  define  the 
expectation  value  on  core  m by  restricting  the  orbitals  fi  to  be  only  those  on  core 
m.  It  can  be  shown,  applying  this  formulation  to  the  orientation  parameters  in  the 
chosen  geometry,  that  they  are  all  identically  equal  to  zero. 


This  would  not  necessarily  be  true  if  the  collision  geometry  were  such  that  the 
collision  occurred  in  the  x-y  plane.  Therefore,  some  of  the  workers  in  the  field 
report  orientations  while,  because  of  our  choice  of  geometry,  we  will  only  report 
an  alignment  parameter.  In  our  geometry  the  only  alignment  parameter  with 
nonzero  value  is 


We  will  be  interested  in  reporting  what  is  called  the  integral  alignment,  A20  for 
comparison  to  other  theoretical  and  experimental  results.  The  integral  alignment 
is  defined  for  a given  projectile  energy  E 


= {Jz)m  = 0 


(5.8) 


V6 


(5.9) 


A2o{E)  = 


(Tx{E)  + az{E) 


(5.10) 


where  ctx{E)  and  (Tz{E)  are  the  state-to-state  integral  cross  sections  at  projectile 
energy  E for  the  2px  and  2pz  atomic  orbitals  at  core  m,  respectively.  The  state- 
to-state  integral  cross  section  to  a final  atomic  orbital  fi  for  projectile  energy  E 
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is  defined  by 


00 


(5.11) 


0 


where  b is  the  impact  parameter  and  P^q{E,  b)  is  the  transition  probability  from 
initial  state  a to  final  state  fx  at  the  final  time. 

Another  concept  from  time-independent  studies  that  will  prove  to  be  of  interest 
for  time-dependent  studies  is  the  idea  of  atomic  population.  It  has  been  shown 
within  a time-independent  framework  that  the  electronic  population  of  an  atomic 
state  may  be  generally  defined  by  [48] 


where  the  choice  of  a is  restricted  to  the  interval  0 < q < 1.  We  will  consider 
how  these  definitions  perform  as  measures  of  time-dependent  populations.  The 
two  particular  choices  of  a with  which  we  will  concern  ourselves  define  the 
so-called  Mulliken  and  LOwdin  population.  The  Mulliken  population  is  defined 


a = i.  The  calculation  of  the  Lowdin  population  involves  finding  the  eigenvalues 
of  the  overlap  matrix  and  hence  requires  a little  more  care  in  calculation  than  the 
Mulliken  population;  however,  the  calculation  of  populations  requires  very  little 
of  the  actual  computational  time  for  any  given  collision  so  that  time  involved  in 
population  calculation  will  not  be  considered  critical  to  the  evaluation  of  differing 
population  analyses. 


(5.12) 


by  the  choice  of  q = 0 while  the  Lowdin  population  is  defined  by  the  choice  of 
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The  physical  implications  of  the  population  analyses  will  be  the  principal 
criterion  for  this  evaluation.  It  can  be  seen  that  we  may  now  define  the  time- 
dependent  alignment  parameter  in  terms  of  the  Mulliken  population 

A^^\t)  = N2p,{t)-2N2p,{t)  (5.13) 

This  definition  of  alignment  is  seen  to  be  an  absolute  rather  than  relative  definition 
with  regard  to  the  total  2p  population.  In  summary,  these  definitions  of  time- 
dependent  properties  have  been  chosen  because  of  their  close  relationship  to 
already  familiar  properties  calculated  in  time-independent  approaches.  It  is  hoped 
that  these  time-dependent  measures  of  atomic  populations  and  electronic  alignment 
will  provide  new  insights  into  the  most  fundamental  electronic  rearrangement 


interactions. 


CHAPTER  6 

INTEGRATION  OF  THE  LINEARIZED  TDHF  EQUATIONS 


So  far  we  have  considered  the  derivation  and  development  of  a temporally 
linearized  time-dependent  Hartree-Fock  description  for  the  electronic  density 
matrices  in  collision  involving  electronic  rearrangement.  In  this  chapter  we  will 
describe  the  calculational  approach  to  the  solution  of  the  first-order  differential 
equation  in  time  previously  derived  for  the  motion  of  the  cores,  which  are  treated 
as  classical  particles  within  the  eikonal  approximation,  and  for  the  electronic 
densities,  which  are  treated  as  quantum  mechanical  entities.  It  will  be  seen  that 
the  order  in  which  the  integrations  are  accomplished  is  essential  to  achieving 
the  appropriate  result.  Further,  the  procedure  for  including  the  coupling  between 
electronic  and  nuclear  motions  will  be  discussed  in  greater  detail. 

In  chapter  4,  we  began  the  discussion  of  the  criteria  for  the  choice  of 
an  appropriate  time  step  developed  with  particular  consideration  given  to  the 
electronic  degrees  of  freedom.  We  recall  that  the  time  step  At  is  chosen  so  that 
the  condition 


Slower  ^ ||p7||  “ ^higher 


(6.1) 


is  satisfied.  The  propagation  of  the  core  degrees  of  freedom  require  the  solution 
of  the  Hamilton  s equations,  which  need  the  knowledge  of  the  core  positions  and 
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momenta  at  a given  time  and  the  potential  in  which  the  cores  move.  We  have 
chosen  to  use  the  Runge-Kutta  fourth  order  algorithm  for  the  integration  of  these 
Hamilton’s  equations.  It  is  important  that  the  Runge-Kutta  fourth  order  algorithm 
is  self  starting  and  can  be  adapted  to  use  our  choice  of  time  step  rather  than  needing 
a constant  time  step.  This  choice  of  integrator  is  by  no  means  unique  and  perhaps 
may  not  be  optimal,  however  for  the  purposes  of  our  calculation  the  algorithm 
is  known  to  be  quite  reliable.  The  integration  of  the  core  coordinates  requires 
a small  fraction  of  the  total  computing  time  involved  in  a typical  computation 
so  that  we  have  not  put  a considerable  amount  of  effort  into  the  optimization  of 
the  core  integrations. 

Once  we  have  chosen  a test  time  step  At  , the  time  step  is  taken  in  two 
parts.  First  the  time  is  advanced  by  one  half  of  the  total  time  step.  As  we 
know  the  full  electronic  density  matrix  at  time  to,  we  may  calculate  the  potential 
which  depends  on  this  density  matrix  at  time  to.  From  the  positions  and  momenta 
of  the  cores  at  time  to  as  well  as  the  potential  due  to  the  electronic  degrees  of 
freedom,  we  are  able  to  move  the  cores  to  their  positions  and  momenta  at  time 
^1/2  = <0  + 5 At.  After  the  core  coordinates  have  been  calculated  at  time  ti/2,  we 
are  able  to  calculate  the  reference  density  at  time  ti/2.  With  this  new  reference 
density,  we  are  now  ready  to  complete  the  time  step.  Completing  the  time  step 
requires  the  calculation  of  the  electronic  potential  at  time  ti/2,  which  is  done  using 
the  new  reference  density.  Next  the  core  coordinates  are  propagated  to  the  end  of 
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the  time  step  ti.  Having  the  core  positions  and  momenta  at  time  ti  we  are  able  to 
use  the  temporal  linearization  scheme  to  calculate  the  reference  density  and  the 
change  in  the  density  at  time  ti.  Schematically  the  procedure  looks  like  this 

Given  (R,P,p,V)/o  — ♦ Advance  to  ti  — > 

2 

Calculate  (R,  P,  ♦ Advance  to  t-[  (6.2) 

Calculate  (R,  P,p,  V)tj 

If  the  condition  on  the  change  of  the  density  matrix  is  not  satisfied  another  time 
step  is  chosen.  If  the  change  of  the  density  matrix  is  too  small  the  original 
time  step  is  tripled  while  if  the  change  of  the  density  matrix  is  too  large  the 
original  time  step  is  halved.  Note  that  these  choices  of  multiples  for  the  time 
step  avoid  the  possibility  of  ’trapping’  in  which  the  change  in  the  density  matrix 
is  alternately  too  large  then  too  small  indefinitely.  If  the  change  in  the  density 
matrix  is  within  tolerances  the  time  step  is  accepted,  time  ti  is  then  renamed 
time  to  and  the  previous  time  step  serves  as  a first  guess  to  the  next  time  step. 
When  the  separation  between  the  cores  exceeds  an  ending  value  the  integration  is 
completed  and  the  final  results  of  the  computation  are  reponed.  The  calculation 
of  a single  integrated  cross  section  may  require  that  the  integration  be  carried 
out  for  as  many  as  100  impact  parameters  and  that  the  transition  probabilities  be 
integrated  over  these  impact  parameters. 

There  is  a considerable  advantage  to  calculating  each  time  step  in  two  halves. 
The  integral  Ti„pj  which  depends  on  the  driving  term  in  equation  4.15  can  be 
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calculated  to  order  At^  by  expanding  the  driving  term  around  time  ti/2  rather 
than  around  time  to  or  time  ti . In  the  calculations  present  here  the  driving  term  is 
expanded  around  time  ti/2  and  terms  of  first  order  in  time  and  of  higher  order  are 
ignored  leaving  only  the  value  of  the  driving  term  itself  at  time  ti/2-  Tests  have 
been  carried  out  by  keeping  the  first  order  term  in  time  as  well,  but  these  tests 
showed  no  considerable  difference  in  the  calculation  when  run  under  the  same 
initial  conditions  and  given  the  same  tolerances. 

A series  of  test  calculations  have  been  run  to  ensure  the  convergence  and 
stability  of  the  numerical  procedures  for  the  integration  of  the  equations  for  core 
and  electronic  density  matrix  propagation.  These  tests  have  involved  varying 
of  some  parameters  input  to  the  calculation  while  holding  others  fixed.  The 
parameters  which  have  been  tested  include: 

• the  initial  and  final  separation  of  the  cores 

• the  higher  and  lower  tolerances 

• the  initial  time  step 

• the  maximum  time  step 

The  choice  of  basis  set  for  the  calculation  also  plays  an  important  role  in  the 
quality  of  the  result,  the  next  chapter  will  be  dedicated  to  the  discussion  of  the 
choice  of  basis  set  and  the  calculation  of  integrals  in  the  various  bases. 


CHAPTER  7 

CHOICE  OF  BASIS  SET  AND  ITS  EFFECT  ON  CALCULATION 


The  expansion  of  the  density,  overlap,  and  Fock  operators  into  matrix  form 
requires  the  choice  of  a set  of  basis  functions  in  whose  terms  the  atomic  and 
molecular  orbitals  may  be  expanded.  The  choice  of  basis  set  has  been  dealt 
with  extensively  in  the  case  of  time-independent  atomic  and  molecular  structure 
calculations.  A time-dependent  dynamics  calculation  presents  a series  of  problems 
that  are  not  necessarily  confronted  in  the  time-independent  case.  The  need  for 
geometries  far  away  from  equilibrium  as  well  as  the  large  range  of  possible 
separations  of  the  cores  is  one  such  problem.  A second  problem  is  the  need  for 
the  inclusion  of  a large  number  of  atomic  and  molecular  states  which  are  not  in 
the  ground  state  configuration.  The  need  in  time-dependent  methods  to  calculate 
overlap  and  Fock  matrix  elements  at  the  large  number  of  differing  geometries 
which  occur  over  the  course  of  the  collision  also  separates  dynamics  and  structural 
calculations.  Finally,  the  need  to  interpret  the  final  result  in  terms  of  the  atomic 
orbitals  of  the  products  of  the  interaction  puts  a premium  on  an  insightful  choice 
of  basis. 


38 


39 


Recall  equation  3.4,  where  we  expanded  the  time-dependent  molecular  orbitals 
as  linear  combinations  of  traveling  atomic  orbitals 


It  is  rather  evident  that  as  the  core  velocity  Vm  tends  to  0 the  value  of  the  ETF 
tends  to  1.  Hence  it  would  be  interesting  to  study  to  effect  of  the  ETF  in  a 
few  simple  cases  to  determine  at  what  core  velocity  it  becomes  appropriate  to 
neglect  ETFs.  In  Figures  7-1  and  7-2  we  display  the  result  of  a calculation  on 
the  proton-Hydrogen  atom  system  using  the  traveling  atomic  basis  and  the  non- 
traveling, static  atomic  basis.  Here  we  have  used  a basis  which  includes  only 
the  Is  hydrogenic  orbital  on  each  center  and  the  screened  Coulomb  potential,  as 
described  in  chapter  2,  is  used  for  the  calculation  of  the  nuclear  trajectories. 

The  transfer  probability  is  plotted  versus  the  impact  parameter  at  impact 
energies  of  1 keV  and  10  eV,  respectively.  It  is  clear  that  the  result  at  1 keV 
which  corresponds  to  a velocity  of  0.2  au  is  profoundly  affected  by  the  use  of 
ETFs.  The  10  eV  result  which  corresponds  to  a velocity  of  0.02  au  shows  less 
of  a difference  between  TAO  and  AO  bases  as  the  two  are  distinguishable  only 
at  very  small  values  of  the  impact  parameter.  We  conclude  therefore,  that  the 


(7.1) 


and  used  the  local  version  of  the  electron  translation  factor 


(7.2) 
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collisions  at  a lab  energy  of  1 keV.  (full  line;  including  ETFs,  dotted  line;  without  ETFs) 
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collisions  at  a lab  energy  of  10  eV.  (full  line:  including  ETFs,  dotted  line:  without  ETFs) 
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use  of  ETFs  for  collision  velocities  of  0.02  au  or  below  is  not  necessary  in  the 
calculation  of  integrals  if  one  wishes  to  obtain  correct  integral  cross  sections.  It 
is  however  necessary  in  the  formalism  to  cancel  spurious  long-range  couplings. 

We  have  used  two  types  of  bases  with  and  without  ETFs  for  the  representation 
of  the  electron  density  matrices  in  our  time -dependent  studies.  The  first  is  the 
exact  hydrogenic  basis  in  which  we  include  the  Is,  2s,  2px,  2py,  and  2pz  orbitals 
on  each  center.  The  overlap  and  Hamiltonian  elements  are  programmed  as  a 
function  of  the  effective  charge  of  the  core  on  which  the  orbital  is  centered.  The 
integration  of  the  overlap  and  Hamiltonian  matrix  elements  involving  two  centers 
was  carried  out  over  elliptical  coordinates  defined  by 


( = + rbj) 

n - rtj) 

^ = ^a  = n 
dv  = ~ T)‘^)didT)d^ 

where  the  ranges  of  integration  are 


(7.3) 


l<^<oo;  -l<7/<l;0<(y?<27r  (7.4) 


Here  we  have  defined  the  separation  between  centers  a and  b to  be  R,  the  electron 
coordinate  referred  to  center  a to  be  raj,  and  the  electron  coordinate  referred  to 
center  b to  be  rbj.  The  azimuthal  angle  is  taken  to  be  (p  and  the  coordinate  system 
is  oriented  so  that  the  z-axis  is  along  the  intercenter  axis  and  the  x-axis  is  in  the 
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plane  of  the  collision.  The  resulting  integrals  in  this  coordinate  system  need  to 
be  rotated  to  the  frame  in  which  the  time  propagation  is  being  performed.  This  is 
accomplished  by  a rotational  similarity  transformation  among  integrals  involving 
p-type  orbitals  as  given  in  Goldstein’s  Classical  Mechanics,  chapter  6 [49].  The 
rotation  is  performed  around  the  y-axis  through  the  angle  between  the  intercenter 
axis  and  the  time  propagation’s  z-axis.  The  result  of  the  integrations  of  hydrogenic 
orbitals  for  overlap  and  Hamiltonian  matrix  elements  is  presented  in  Appendix  A 
as  a function  of  intercenter  separation  and  effective  core  charge. 

A popular  choice  of  basis  set  for  time-independent  calculations  is  the  basis  of 

Gaussian  type  orbitals.  It  has  been  shown  by  Boys  [50]  that  the  use  of  Gaussian 

type  orbitals  allows  for  the  rapid  calculation  of  one-  and  two-electron  integrals  in 

the  context  of  time-independent  Hartree-Fock  theory.  It  seems  natural  to  employ 

a similar  strategy  in  time-dependent  calculations.  This  strategy  entails  writing  of 

atomic  orbitals  in  the  form  (after  scaling  to  eliminate  charge  parameters)  [51] 

K 

(7.5) 

k 


where  rj  = x,  y,  z.  Here  the  Gaussian-type  orbital  for  52pyand  gop,  are 


, . , ,128q5 

92py\<^i‘>'mj)  — \ ^3 

^ ^_,128a5 

92pz[Ot,  Cmj)  — ( 3 

7T 


)ix„,jexp(-Qr^y) 
)«y^jexp(-Qr^j) 
)hmj  exp(-arij) 


(7.6) 
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where  the  number  K ranges  in  our  studies  from  3 to  6,  are  the  Gaussian 
exponential  coefficients,  and  are  the  expansion  coefficients.  The  Slater-type 
orbitals  are  identical  to  hydrogenic  orbitals  for  the  Is  and  2p  orbitals.  For  the  2s 
orbital  the  use  of  Slater-type  orbitals  written  as  linear  combinations  of  Gaussian- 
type  orbitals  is  not  sufficient  as  the  2s  Slater-type  orbital  is  nodeless.  Therefore, 
Gaussian  exponential  coefficients  and  expansion  coefficients  have  been  chosen  to 
fit  the  hydrogenic  2s  function  rather  than  the  Slater-type  function. 

There  are  a number  of  available  integral  packages  for  the  integration  of 
Gaussian-type  orbitals.  We  have  chosen  to  use  the  MOLECULE  package  devel- 
oped by  Jan  Almloff  and  P.  R.  Taylor  [52].  The  local  VMOL  package  calculates 
both  one-  and  two-electron  integrals  for  either  segmented  contracted  or  generally 
contracted  orbitals.  We  have  chosen  to  use  segmented  contracted  sets  of  Gaussians 
to  represent  the  Is,  2s  and  2p  hydrogenic  orbitals  for  a given  effective  charge. 

The  VMOL  package  writes  the  results  for  a given  geometry  to  integral 
files.  This  would  necessitate  numerous  input/output  operations  in  our  time- 
dependent  code  which  could  be  quite  time  consuming.  Hence,  we  have  modified 
the  MOLECULE  package  to  write  the  results  of  integral  calculations  into  the 
appropriate  arrays  and  made  the  subroutine  calls  internal  to  our  main  program 
so  that  we  are  able  to  calculate  the  integrals  at  the  large  number  of  geometries 
we  require  and  save  the  computing  time  that  would  otherwise  be  devoted  to 
MOLECULE  input/output. 
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We  have  also  tested  the  effect  of  varying  the  number  of  Gaussians  used 
to  approximate  the  hydrogenic  orbitals.  In  Figure  7—3  we  show  the  transfer 
probability  versus  impact  parameter  for  a proton  colliding  with  a hydrogen  atom 
in  its  ground  state. 

This  calculation  was  carried  out  at  an  impact  energy  of  1 keV  using  from  3 to 
6 Gaussians  to  represent  the  Is  orbital  of  hydrogen.  Only  the  Is  orbitals  on  each 
center  were  used  for  this  calculation.  Also  included  is  the  result  using  the  exact 
Is  hydrogenic  on  each  center.  We  see  that  the  agreement  in  general  is  rather  good 
while  there  is  a distinct  improvement  in  moving  from  3 to  4 Gaussians  per  orbital 
the  improvement  beyond  4 is  only  marginal  and  integrated  cross  sections  would 
not  be  markedly  improved.  We  move  now  to  the  consideration  of  inclusion  of 
electron  translation  factors  in  the  one-electron  integrals. 

The  inclusion  of  ETFs  in  two-center,  one-electron  integrals  done  over  hy- 
drogenic orbitals  is  problematic  because  the  integrations  cannot  be  carried  out 
analytically.  Various  numerical  procedures  have  been  implemented  which  require 
the  numerical  integration  over  one  or  two  variables.  Another  option  that  has  been 
tried  in  the  literature  is  to  write  the  orbitals  as  linear  combinations  of  Gaussians. 
In  this  scheme  the  analytical  forms  of  the  integrals  can  be  obtained  and  the  ap- 
proximation occurs  only  in  the  fitting  of  the  Gaussians  to  the  actual  orbitals  of 
interest.  We  have  expanded  the  work  of  McWeeny  by  including  higher  angular 
momentum  states  and  adding  a version  of  the  2s  orbital  similar  in  form  to  the  Is 


1.00 
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Figure  7—3:  Electron  transfer  probability  vs.  Impact  parameter  in  au  (plus:  exact  Is  hydrogenic, 
triangle:  3 Gaussians;  diamond:  4 Gaussians;  box:  5 Gaussians;  cross:  6 Gaussians) 
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orbital  rather  than  McWeeny 


’s  choice  of  2s  orbital  which  is  of  the  form  [53] 


9s{a,rmj)  = (— )«r2  exp(-ar2  ) 


(7.7) 


The  derivation  of  the  higher  angular  momentum  states  is  greatly  aided  by 
the  observation  that  the  Gaussian  basis  functions  of  one  angular  momentum  are 
related  to  the  derivative  of  those  of  lower  angular  momentum  states.  Given  that 

9s{o!,fmj)  = (— )<  exp(-ar^  ) (7.8) 

7T 


if  we  take  the  derivative  of  the  s- state  Gaussian  with  respect  to  a coordinate  of 
the  center  m we  obtain 

-^9s{a,rmj)  = { — )*ar)mexp{-arl,.) 

^ (7.9) 


where  the  index  t)  runs  over  the  x,  y,  and  z coordinates  of  center  m.  Similarly, 
we  have  that  taking  the  derivative  of  the  p-type  Gaussian  with  respect  to  the 
coordinates  of  the  center  m 


-7r-9pr,  ^mj  ) = V^9dr,(  («,  ^mj  ) (7. 1 0) 

It  should  be  mentioned  at  this  point  that  our  time-dependent  program  was  altered 
to  accommodate  complex  orbitals  so  that  the  ETFs  could  be  incorporated.  There  is 
now  an  obvious  scheme  for  the  calculation  of  the  two-center  overlap  integrals  in 
which  they  are  related  to  s-type  overlaps  with  differing  Gaussian  exponential 
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coefficients  and  coefficients  derived  by  differentiation.  The  results  of  these 
operations  appear  in  Appendix  B. 

We  now  need  to  turn  our  attention  to  the  problem  of  calculating  two-center 
Hamiltonian  matrix  elements  including  RTFs.  The  most  important  difficulty,  that 
of  calculating  the  core-electron  interaction,  has  been  addressed  by  other  authors. 
Their  approaches  rely  on  a number  of  approximations  that  we  either  are  unable  to 
make  or  would  prefer  not  to  make.  We  will  avoid  the  problem  of  the  calculation 
of  the  core-electron  interaction  by  making  the  following  observation.  We  write 
the  Hamiltonian  operator  in  the  following  suggestive  form 

H = f + Vm  + t-\-Vr.-t  (7.11) 

where  we  have  added  and  subtracted  the  kinetic  energy  operator  T.  Taking  the 
matrix  element  of  H between  traveling  atomic  orbitals  on  two  differing  centers 
we  obtain 


= (7.12) 


where  tfi  is  the  eigenenergy  of  the  orbital  for  its  appropriate  asymptotic 
Hamiltonian 


H^  = T + Vm 


(7.13) 


The  problem  of  calculating  two-center,  one-electron  Hamiltonian  elements  now 
reduces  to  the  more  tractable  problem  of  calculating  the  kinetic  energy  matrix 
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elements  which  may  in  turn  be  related  to  the  two-center,  one-electron  overlap 
integral  previously  discussed.  The  results  of  these  derivations  are  also  presented 
in  Appendix  B. 

Testing  of  the  two-center  overlap  and  one-electron  Hamiltonian  elements  has 
included  the  comparison  of  traveling  atomic  orbital  calculation  to  static  atomic 
orbital  calculation  in  the  zero  velocity  limit.  Agreement  to  at  least  9 decimal 
places  was  found  between  our  calculations  and  calculations  made  for  the  same 
Gaussian  basis  using  the  MOLECULE  program  package.  Further  testing  involved 
the  calculation  of  overlap  and  Hamiltonian  matrix  elements  for  a small  basis  and 
checking  these  results  by  hand.  These  results  agreed  to  at  least  7 decimal  places, 
where  the  hand  calculations  were  restricted  by  the  accuracy  of  the  calculator  used. 


CHAPTER  8 

APPLICATIONS  I:  THE  PROTON-HYDROGEN  ATOM  SYSTEM 


Overview  of  Experimental  Methods 


We  first  applied  our  method  to  the  simplest  of  collisional  systems  consisting 
of  two  protons  and  one  electron,  that  is  the  collision  of  an  incident  proton  on  a 
target  Hydrogen  atom.  Despite  its  apparent  simplicity,  this  system  has  not  been 
solved  exactly  and  hence,  has  inspired  a great  deal  of  theoretical  and  experimental 
study.  Properties  of  interest  for  these  collisions  are  the  total  integral  cross 
section,  the  state-to-state  integral  cross  section,  and  the  differential  cross  section 
for  electron  transfer,  electron  transfer  excitation,  and  direct  electron  excitation. 
Integral  alignment,  as  described  in  chapter  5,  has  also  been  of  theoretical  and 
experimental  interest.  Energies  of  interest  for  the  incident  proton  have  ranged 
from  as  low  as  1 eV  to  greater  than  1 MeV.  We  concern  ourselves  with  the  lower 
end  of  this  range,  that  is  2 keV  and  below.  We  will  discuss  the  role  of  electron 
translation  factors,  basis  set  type  and  size,  and  choice  of  electronic  potential  in 
regard  to  their  effect  on  the  calculation  of  the  above  quantities  of  interest.  We 
will  also  discuss  the  relative  difficulties  of  correct  calculation  of  these  quantities 
over  the  entire  energy  range  of  interest. 
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The  numerous  experiments  that  have  been  carried  out  on  this  system  can  be 
divided  into  five  categories  which  depend  on  the  way  the  protons  interact  with 
the  Hydrogen  atoms  and  the  way  the  final  states  are  detected.  There  are  two 
experimental  geometries  for  the  interaction 

• The  high  energy  proton  beam  is  incident  on  a thermal  energy  target  Hydrogen 
atom  beam,  a crossed  beam  geometry. 

• The  proton  beam  and  the  target  Hydrogen  atom  beam  travel  through  an 
interaction  region  together,  a merging  beam  geometry. 

Detection  apparati  vary  according  to  the  quantity  to  be  measured.  There  are  four 
detector  arrays  considered  here 

• An  ion  or  atom  counting  mechanism  collects  the  forward  scattered  projectiles 
and  sorts  them  according  to  their  charge,  a counter. 

• A polarization  sensitive  photon  counter  collects  radiation  given  off  by  either 
projectile  or  target  atoms,  a Lyman-o  detector. 

• A quenching  electric  field  is  applied  to  the  collision  region  and  radiation 
collected,  a quenched  Lyman-o  detector. 

• An  ion  counter  and  an  atom  counter  are  used  to  analyze  projectile  beams 
selected  by  their  scattering  angle,  a differential  detector. 

The  combination  of  a crossed  beam  geometry  and  a counter  detector  allows 
the  measurement  of  total  integral  cross  sections  at  fairly  high  energies  [54-57]. 
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The  combination  of  a crossed  beam  geometry  and  a Lyman-a  detector  allows 
the  measurement  of  state-to-state  integral  cross  section  to  2p  states  and  integral 
alignments  [58,  46,  59,  60].  The  combination  of  a crossed  beam  geometry  and 
a quenched  Lyman-a  detector  allows  the  measurement  of  state-to-state  integral 
cross  sections  to  2p  and  2s  states  [61-64].  The  combination  of  a crossed  beam 
geometry  and  a differential  detector  allows  the  measurement  of  differential  cross 
sections  [65].  Finally,  the  combination  of  a merging  beam  geometry  and  a counter 
detector  allows  the  measurement  of  total  integral  cross  sections  at  lower  energies 
[66]. 


Integral  Cross  Sections 


The  determination  of  total  integral  cross  section  at  higher  energies  is  a decep- 
tively easily  calculated  quantity.  A comparison  of  theoretical  and  experimental 
values  obtained  in  the  case  of  a 1 or  2 keV  projectile  incident  on  Hydrogen  atom 
is  shown  in  table  8-1.  Table  8-1  illustrates  that  the  total  integral  cross  section  at 
an  intermediate  energy  is  quite  easily  calculated  to  an  accuracy  of  20%  or  less. 
This  raises  the  question  of  why  this  should  be  true  even  for  theoretical  methods 
which  use  only  the  Is  orbital  on  each  center.  A more  detailed  study  of  the  role  of 
basis  set  size  and  type  and  the  role  of  the  potential  in  which  the  nuclei  move  may 
enhance  our  insight  as  to  why  this  should  be.  Figure  8-1  shows  the  influence  of 
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the  choice  of  nuclear  trajectory  on  a calculation  at  1 keV  with  a Is  hydrogenic 
orbital  per  center  . 

It  is  obvious  that  the  choice  of  nuclear  trajectory  does  not  impact  the  cal- 
culational  result  in  the  case  of  the  1 keV  collision.  The  question  becomes:  is 
the  choice  of  nuclear  trajectory  ever  important  for  the  calculation  of  total  inte- 
gral cross  sections?  It  becomes  apparent  that  the  choice  of  nuclear  trajectory  is 
likely  to  be  more  important  as  the  projectile  energy  is  lowered.  The  effect  of 
the  screened  Coulomb  potential  is  expected  to  allow  for  some  curvature  in  the 
trajectory  but  not  as  much  curvature  as  would  be  expected  from  the  Coulomb 
potential.  We  therefore  present  in  Figure  8-2  a comparison  of  nuclear  trajectories 
at  an  impact  energy  of  100  eV.  Again,  the  calculation  includes  a Is  hydrogenic 
orbital  per  center. 

Here  we  see  that  at  the  lower  energy  of  100  eV  the  zero  potential  is  closer  to 
the  screened  Coulomb  potential  than  the  Coulomb  potential.  It  is  anticipated  that 
the  screened  Coulomb  potential  will  do  the  best  job  in  calculating  total  integral 
cross  sections  at  lower  energies.  We  will  come  back  to  this  assertion  in  good 
order.  We  now  continue  to  examine  the  role  of  the  choice  of  nuclear  trajectory 
at  lower  energies  by  comparing  these  choices  at  10  eV.  In  Figure  8-3  we  show  a 
comparison  of  the  results  for  these  trajectory  choices  at  an  impact  energy  of  10  eV. 

Here  we  see  that  the  three  trajectories  choices  give  very  different  results  with 
the  Coulomb  potential  departing  most  from  the  screened  Coulomb  potential.  A 


Table  8.1:  A comparison  of  experiment  and  theory  for  total  integral  cross  sections. 
Total  Integral  Cross  Section  x 10*^*  cm^  Theory  and  Experiment 
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stricter  test  of  the  quality  of  the  trajectory  is  the  comparison  with  differential  cross 
sections  measured  by  experiment  which  we  will  deal  with  later  in  this  chapter. 
Now  let  us  complete  our  discussion  of  total  integral  cross  section,  in  Figure  8-4 
we  show  a comparison  of  total  integral  cross  section  between  our  results  and 
experiments. 

Note  that  the  agreement  between  calculated  and  experimental  results  is  good 
through  the  entire  range  of  projectile  energies  from  2 to  2000  eV.  This  is  another 
bit  of  evidence  that  the  screened  Coulomb  potential  is  the  best  choice  for  this 
system  as  a poorer  choice  of  potential  would  not  give  good  agreement  over  so 
wide  a range  of  energies,  particularly  at  lower  energies. 

We  now  move  on  to  state-to-state  integral  cross  sections.  The  role  of  electron 
translation  factors  in  these  calculations,  which  was  discussed  in  chapter  7,  is  also 
important  in  obtaining  accurate  values  for  state-to-state  integral  cross  sections. 
ETFs  are  most  important  in  the  calculation  of  integral  alignments  which  will  be 
discussed  below.  Table  8-2  shows  a detailed  comparison  of  various  experimental 
and  theoretical  results  in  a 2 keV  impact  of  a proton  on  a target  Hydrogen  atom. 
This  comparison  shows  that  our  theoretical  results  are  in  good  agreement  with 
experiment  for  all  state-to-state  integral  cross  sections.  Note  that  the  predominate 
mechanism  for  charge  transfer  is  the  resonant  transfer  channel,  this  helps  to  explain 
why  total  electron  transfer  integral  cross  sections  are  so  easily  found  to  within  20% 
theoretically.  The  most  difficult  state-to-state  cross  sections  to  obtain  accurately 
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Table  8^:  Comparison  of  theoretical  and  experimental  results  for  state-to-state  integral  cross  sections  for  a 2 keV  proton  incident  on  a Hydrogen  atom. 
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are  the  2s,  2px,  and  2pz  which  require  the  inclusion  of  ETFs  in  our  method.  Note 
that  all  other  theoretical  methods  in  this  comparison  include  ETFs  in  one  way  or 
another.  Another  experimentally  measurable  quantity  that  allows  us  to  test  the 
sensitivity  of  our  results  is  the  integral  alignment.  Accurate  calculation  of  integral 
alignments  requires  a basis  set  which  includes  ETFs  and  has  sufficient  size  to  give 
a good  description  of  the  n=2  orbitals  of  the  Hydrogen  atom.  In  this  study  of 
integral  alignment  our  basis  set  has  consisted  of  5 hydrogenic  orbitals  per  center 
each  written  as  a sum  of  6 Gaussian  orbitals.  The  5 hydrogenic  orbitals  that 
were  used  include  the  Is,  2s,  2px,  2py,  and  2pz.  Table  8-3  shows  a comparison 
of  theoretical  and  experimental  results  for  integral  alignments  for  a number  of 
energies  in  the  proton- Hydrogen  atom  system. 

Here  the  result  of  Fite  et  al.  [59]  appears  to  be  anomalous.  Our  result  agrees 
well  with  Hippler’s  [46]  at  both  available  energies.  We  also  predict  a drop  in 
integral  alignments  at  lower  energies,  a prediction  made  by  Deumens  et  al.  [39] 
as  well,  however  the  exact  energy  at  which  the  drop  in  integral  alignment  occurs 
is  not  agreed  upon,  Deumens  et  al.  [39]  see  a drop  below  500  eV.  It  is  appropriate 
to  keep  in  mind  that  at  lower  energies  the  excitation  cross  section  amount  to  less 
than  0.1%  of  the  resonant  transfer  cross  section  so  that  the  integral  alignments  at 
these  lower  energies  will  be  very  difficult  to  measure  experimentally. 
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Differential  Cross  Sections 

The  final  experimental  verification  of  the  soundness  of  our  theoretical  method 
comes  from  the  examination  of  differential  cross  sections.  A comparison  of  the 
effect  of  choice  of  potential  on  differential  cross  section  is  shown  in  Figures  8-5 
and  8-6.  We  introduce  the  reduced  differential  cross  section  p as  defined  by 
Houver  et  al.  [65] 


p = a9is'm9L  (8.1) 

where  a is  the  differential  cross  section  in  au  and  9^  is  the  lab  scattering  angle 
in  degrees.  Figure  8-5  shows  the  elastic  reduced  differential  cross  section  versus 
lab  angle  in  degrees  for  the  screened  Coulomb  potential,  the  Coulomb  potential 
and  the  the  experimental  result  of  Houver  et  al.  [65]  at  a projectile  energy  of 
410  eV.  Figure  8-6  shows  the  charge  transfer  reduced  differential  cross  section 
versus  lab  angle  in  degrees. 

It  is  necessary  to  note  that  the  experiment  by  Houver  et  al.  [65]  has  an  angular 
resolution  of  0.2  degrees  at  small  angles  but  at  larger  angles  the  resolution  drops 
to  1.2  degrees.  Hence  it  is  quite  likely  that  the  calculated  minimum  in  the  elastic 
reduced  differential  cross  section  at  about  4 degrees  would  be  missed  by  the 
experiment.  There  is  substantial  agreement  between  our  work  and  the  experiment 
in  term  of  the  maxima  around  1.5  degrees  and  3.0  degrees  and  the  large  angle  trend 
around  7 degrees  seems  to  show  good  agreement.  The  calculation  of  differential 
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incident  on  a Hydrogen  atom  target  {full  line:  SCP;  dotted  line:  CP;  triangles:  [65]  } 
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cross  section  vs.  scattering  angle  is  quite  sensitive  to  the  choice  of  effective 
potential.  Choosing  the  Coulomb  potential  leads  to  scattering  angles  which  are 
about  twice  as  large  as  those  calculated  using  the  screened  Coulomb  potential. 
Obviously,  the  use  of  the  zero  potential  is  meaningless  as  it  leads  asymptotically 
to  a zero  scattering  angle  for  any  impact  parameter.  We  conclude,  therefore  that 
the  choice  of  the  screened  Coulomb  potential  is  essential  for  accurate  calculation 
of  differential  cross  section  at  all  energies  of  interest.  We  continue  to  compare  our 
results  using  the  screened  Coulomb  potential  with  those  of  experiment.  It  should 
be  noted  that  Houver  used  a “best  fit”  to  other  theoretical  results.  It  is  not  suprising 
that  our  results  show  a consistent  0.4  degree  difference  in  the  positions  of  extrema 
with  their  results.  The  most  important  features  for  comparison  are  the  general 
features  of  the  reduced  differentail  cross  sections  and  the  difference  in  angle 
between  maxima  and  minima.  Figures  8-7  and  8-8  present  a comparison  of  our 
calculated  results  using  the  screened  Coulomb  potential  and  those  of  experiment 
for  elastic  and  charge  transfer  reduced  differential  cross  sections  versus  lab  angle 
in  degrees  for  500  eV  and  700  eV  projectiles,  respectively. 

Again  we  note  the  agreement  with  the  experimental  results  with  respect  to 
the  type  and  spacing  of  the  extrema  of  the  reduced  differential  cross  sections. 
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Time-dependent  Properties 


We  now  move  on  to  consider  the  time  dependent  measures  of  atomic  pop- 
ulation discussed  in  chapter  6.  We  begin  with  a comparison  of  Mulliken  and 
Ldwdin  populations.  Figures  8-9  and  8-10  present  a comparison  of  the  time 
evolution  of  the  Mulliken  and  Ldwdin  population  in  a 1 keV  proton-Hydrogen 
atom  collision,  respectively.  The  impact  parameter  is  0.2  au,  the  zero  of  time  is 
chosen  to  occur  at  the  distance  of  closest  approach,  and  the  screened  Coulomb 
potential  has  been  used. 

Note  that  the  Mulliken  population  is  occasionally  greater  than  1 or  less  than 
0.  This  is  not  a behavior  we  desire  in  our  measure  of  atomic  population.  As  the 
Ldwdin  population  does  not  show  these  pathologies  we  will  prefer  the  Ldwdin 
population  as  a measure  of  atomic  population  as  a function  of  time.  Figures  8-11 
through  8-17  show  the  time  evolution  of  the  Ldwdin  populations  of  the  other 
target  and  projectile  states  involved  in  this  calculation.  Note  that  the  2py  orbital 
is  not  shown  as  it  is  never  populated  by  the  chosen  geometry  of  the  calculation. 

Notice  that  all  the  populations  oscillate  before  finally  coming  to  a final 
asymptotic  value  which  depends,  of  course,  on  the  energy  of  the  collision  and 
the  impact  parameter.  This  is  an  example  of  the  type  of  physical  insight  that 
we  wish  to  obtain  from  time-dependent  studies.  This  insight  is  again  illustrated 
in  Figures  8-18  and  8-19  in  which  the  target  Is  population’s  time  evolution  is 


1.20 


70 


3 

a 

CA 

d 


8 

8 


CO 


VO 

00 

<s 


o 

<N 


<N 


VO 

00 

ri 

t 


CO 


8 

8 


o 

00 

s 

o 

o 

cs 

8 

d 

d 

d 

d 

9 

9 

Figure  8-9:  Mulliken  population  of  the  target  Is  state  vs  time  in  au  in  a 1 keV  proton-Hydrogen  atom  collision,  at  b=( 
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shown  for  a 10  eV  proton-Hydrogen  atom  collision  for  impact  parameters  of  1.0 
au  and  1.2  au,  respectively. 

Here  we  see  multiple  oscillations  in  the  target  atom  population  as  a function 
of  time.  We  also  notice  that  a slight  change  in  the  impact  parameter  causes  a 
large  change  in  the  result  of  the  collision,  in  this  case  a change  of  0.2  au  in  the 
impact  parameter  makes  the  difference  between  complete  retention  of  the  electron 
by  the  target  and  complete  transfer  to  the  projectile.  Because  of  the  low  collision 
energy  there  is  no  direct  or  transfer  excitation. 

The  final  time-dependent  quantity  that  we  investigate  is  the  alignment  pa- 
rameter derived  in  chapter  6.  To  study  the  time  evolution  of  the  alignment  of 
the  target  we  need  to  consider  a fairly  high  collision  energy.  In  Figure  8-20  we 
present  the  time  evolution  of  the  alignment  in  a 1 keV  collision  at  an  impact 
parameter  of  0.2  au. 

Here  we  see  that  the  alignment  parameter  also  oscillates  during  the  collision 
and  is  actually  quite  negative  in  the  midst  of  the  collision  despite  the  fact  that  the 
alignment  parameter  for  this  collision  is  positive. 

In  summary,  we  have  seen  in  this  chapter  that  the  total  integral  cross  section 
at  fairly  high  energies  is  the  least  demanding  quantity  to  be  accurately  calculated. 
The  differential  cross  section  requires  a realistic  choice  of  electronic  potential 
for  accurate  determination.  The  state-to-state  integral  cross  section  requires 
a traveling  atomic  basis  for  accurate  calculation.  Finally,  we  have  seen  that 
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Figure  8-18:  Lowdin  population  of  the  target  Is  state  vs  time  in  au  in  a 10  eV  proton-Hydrogen  atom  collision,  at  b=1.0  au 
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physical  insight  can  be  gained  by  examining  atomic  populations,  where  the 
Lbwdin  population  is  preferred  over  the  Mulliken  population,  and  by  examining 
the  alignment  parameter  vs.  time. 
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Figure  8-20:  Alignment  parameter  vs  time  in  au  in  a 1 keV  proton-Hydrogen  atom  collision,  at  b=0,2  au 


CHAPTER  9 

APPLICATIONS  II:  THE  HELIUM(2+)-HYDROGEN  ATOM  SYSTEM 


A second  system  of  interest  for  us  is  the  collision  of  a fully  stripped  helium 
ion  with  a target  Hydrogen  atom.  This  system  is  an  example  of  near  resonant 
collisions  in  which  the  electron  transfer  is  predominately  from  the  H Is  state  to 
the  He"^  2s  or  2p  state  which  are  of  equivalent  asymptotic  energy  as  the  original 
system.  As  the  nuclei  are  not  of  the  same  type  the  collision  is  not  truly  resonant 
as  is  the  case  considered  in  chapter  8.  We  will  examine  the  projectile  energy 
range  from  1 keV  to  8 keV.  Again  the  integral  charge  transfer  cross  section  and 
state-to-state  integral  cross  sections  will  be  of  interest  as  well  as  time-dependent 
properties. 

Experiments  on  this  system  can  be  described  in  terms  of  the  categories 
introduced  in  chapter  8.  The  combination  of  crossed  beam  geometry  and  counter 
detection  allows  the  measurement  of  total  charge  transfer  integral  cross  section 
[67].  The  combination  of  crossed  beam  geometry  and  counter  and  Lyman-a 
detection  allows  the  measurement  of  total  charge  transfer  integral  cross  section, 
direct  measurement  of  2p  integral  cross  section  and  indirect  measurement  of  2s 
integral  cross  section  [68-70].  The  combination  of  crossed  beam  geometry  and 
counter,  Lyman-o  and  quenched  Lyman-a  detection  allows  the  measurement  of 
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total  charge  transfer  integral  cross  section  and  direct  measurement  of  2p  and  2s 
integral  cross  sections  [71]. 

We  have  used  a basis  which  includes  the  Is,  2s,  2px  and  2pz  hydrogenic 
orbitals  on  both  the  projectile  and  target  nuclei.  All  calculations  have  been  made 
using  the  zero  potential,  as  the  energy  range  is  high  enough  that  the  calculation 
with  the  average  effective  potential  would  not  improve  the  final  integral  cross 
sections  significantly.  ETFs  have  been  included  by  writing  the  hydrogenic  orbitals 
as  linear  combinations  of  6 Gaussian  orbitals  each.  We  have  investigated  the 
projectile  energies  of  1 keV,  3 keV,  and  8 keV.  Table  9-1  presents  a comparison 
of  our  results  with  various  experimental  and  theoretical  results  at  these  energies. 

The  agreement  with  experiment  is  quite  good  although  our  results  are  consis- 
tently lower  than  some  of  the  other  accurate  theories.  The  lack  of  the  inclusion  of 
the  n=3  shell  orbitals  in  our  calculations  may  help  explain  this  discrepancy.  The 
results  due  to  Rapp  [10]  are  rather  puzzling  as  he  uses  hydrogenic  orbitals  which 
include  ETFs  and  his  approach  gave  good  agreement  with  ours  in  the  proton- 
Hydrogen  atom  system.  The  calculation  of  state-to-state  integral  cross  sections 
is  a more  demanding  exercise.  Table  9-2  presents  a comparison  of  theory  and 
experiment  for  state-to-state  integral  cross  sections  at  a projectile  energy  of  8 keV. 

The  agreement  between  our  state-to-state  integral  cross  section  results  is  again 
much  better  than  those  of  Rapp  [10].  The  slight  difference  between  our  results 
and  other  accurate  theories  may  also  be  a result  of  the  differing  basis  sizes,  as 
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above.  We  see  that  the  agreement  between  theoretical  results  and  experiment  are 
really  quite  good.  Comparison  of  state-to-state  integral  cross  section  at  the  lower 
energies  is  not  presented  due  to  a lack  of  experimental  data. 

Time-dependent  phenomena  that  are  of  interest  for  this  system  include  the 
alignment  of  the  projectile  and  the  electronic  population  of  the  target  Is  state  and 
the  projectile  2s,  2px  and  2p^  states.  Figures  9-1  through  9-4  show  the  Ldwdin 
populations  of  the  target  and  projectile  states.  The  conditions  in  this  calculation 
have  been  chosen  to  correspond  to  the  results  shown  in  the  previous  chapter  for 
the  proton-Hydrogen  atom  system.  This  system,  therefore,  is  the  collision  of  a 4 
keV  He^"^  with  a target  hydrogen  at  an  impact  parameter  of  0.2  au.  Notice  that 
at  this  energy  and  impact  parameter  the  charge  exchange  occurs  primarily  to  the 
2pz  state.  By  comparison  to  the  previous  chapter’s  results  we  see  that  there  is 
less  oscillation  of  the  electronic  population  in  the  current  near  resonant  collision 
than  there  was  in  the  previous  resonant  case.  However,  there  still  is  significant 
oscillation  particularly  in  the  Is  target  and  2pz  projectile  states.  The  variation  in 
time  of  the  alignment  parameter  in  shown  in  Figure  9-5.  This  parameter  shows 
pronounced  oscillations  much  like  those  shown  in  the  previous  chapter.  In  this 
case,  the  alignment  parameter  is  more  negative  than  in  the  proton-Hydrogen  atom 
collision.  Again,  we  gain  considerable  insight  into  the  charge  transfer  mechanism 
as  a function  of  time  in  the  collision  of  an  ion  with  an  atom. 
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CHAPTER  10 

APPLICATIONS  m:  THE  HYDROGEN  ATOM-HYDROGEN  ATOM  SYSTEM 


A third  system  of  interest  for  us  is  the  collision  of  two  Hydrogen  atoms  in  the 
ground  state.  This  is  the  first  two-electron  system  which  we  have  investigated. 
The  total  excitation  integral  cross  section  and  the  state-to-state  integral  cross 
section  will  be  of  interest  as  well  as  time-dependent  properties.  We  will  consider 
the  projectile  energy  of  1 keV. 

Experimental  approaches  to  this  system  can  again  be  described  in  the  terms  of 
chapter  8.  The  combination  of  a crossed  beam  geometry  and  a quenched  Lyman-a 
detector  allows  for  the  measurement  of  the  2s  excitation  and  2s  excitation  transfer 
integral  cross  sections  [63,  61].  The  experimental  results  shown  are  obtained  from 
linear  extrapolation  of  results  obtained  at  2 keV  and  3 keV  projectile  energies. 

In  the  case  of  two  active  electrons  the  electron  spin  becomes  important  and 
the  density  matrices  for  alpha  spin  and  beta  spin  must  be  propagated  through 
time.  So  equations  4.6  and  4.8  become 


and 


zPo"(0=  W"(<o)P^(f)  - P?(f)W“(fo)^ 
tPl(t)=W%)P^(t)-PS(t)W%)^ 

iQ^(t)=  W“(<o)Q"(0  - Q“(f)W“(to)4D“(0 
W^(fo)Q^(0  - Q^(OW^(/o)4d^(0 


(10.1) 


(10.2) 


94 


95 


where  „ 

W"  = S-^F"[P?,Pj] 


(10.3) 


W/J  = S-^F^[P?,Pj] 

A number  of  possible  reactions  must  be  considered  in  the  calculation  of  the 

excitation  cross  sections  in  the  Hydrogen  atom-Hydrogen  atom  system.  These 

possibilities  reflect  the  lack  of  selection  of  the  initial  spin  state  and  of  analysis 

of  the  final  spin  state  in  experimental  results,  and  the  fact  that  the  projectile  and 

target  nuclei  are  identical.  The  electronic  spin  states  could  proceed  from  an  initial 

state  in  which  the  projectile  has  an  electron  in  either  the  alpha  spin  state  or  the 

beta  spin  state  to  a final  state  in  which  the  projectile  has  an  electron  of  either 

the  same  spin  state  or  opposite  spin  state  as  the  initial  state.  The  possible  spin 

states  of  the  initial  states  and  the  final  states  to  the  right  of  the  arrow  may  be 

summarized  by  the  transitions 

aa  -{■  ba  —*  aa  + ba 


-\-  bd  al5  6/3 


aa  4-  6/3  — ♦ <za  -|-  6yS 
—*  a/3  + ba 
al3  + ba  a/3  + ba 


(10.4) 


—+00  + 6/3 

where  a represents  the  projectile  nucleus  and  b represents  the  target  nucleus. 
Therefore,  to  obtain  the  appropriate  excitation  cross  section  for  target  hydrogen 
atom  we  must  average  cross  sections  over  initial  spin  states  and  add  over  final 
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spin  states.  In  particular,  for  electronic  excitation  we  must  add  the  cross  section 
a{aa  + b/3  ^ aa  + b/3),  for  excitation  to  the  final  state  with  the  same  electron 
spins  as  the  initial  spin  state,  to  the  cross  section  cr(ao  + /j/?  a/3  + ba),  for 

excitation  to  the  final  state  with  opposite  electron  spins  at  the  nuclei. 

Table  10-1  shows  a comparison  of  experimental  and  theoretical  results  for 
excitation  cross  sections  to  2s  and  2p  states  of  the  target  Hydrogen  atom  in  the 
collision  of  a 1 keV  projectile  Hydrogen  atom  and  a target  Hydrogen  atom  each 
in  their  ground  electronic  state.  We  use  a basis  including  Is,  2s,  2px,  and  2pz 
Hydrogenic  orbitals  are  written  as  linear  combinations  of  6 Gaussian  orbitals  for 
each  center.  The  electron  translation  factors  are  not  included  in  the  calculation 
of  electron  integrals.  However,  due  to  the  cancellation  of  the  so-called  spurious 
coupling  by  the  ETF’s  in  the  formalism,  the  electronic  populations  do  go  to 
constant  values  a the  final  time.  The  straight  line  trajectory  is  used  insofar  as 
previous  work  shows  that  the  calculation  of  integral  cross  sections  would  not  be 
affected  by  the  use  of  a zero  potential  instead  of  the  averaged  effective  potential 
at  this  energy. 

We  see  in  Table  10-1  that  our  2s  state-to-state  excitation  integral  cross  section 
result  agrees  well  with  the  experiments  and  the  theoretical  results  of  Shingal  et 
al.  [12]  The  disagreement  with  the  result  of  McLaughlin  and  Bell  [11]  is  quite 
sizable  and  raises  questions  regarding  their  procedure.  In  general,  it  is  fair  to  say 
that  the  agreement  of  our  result  with  those  of  other  theories  is  reasonable.  At 
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the  same  time,  the  agreement  with  experiment  seems  quite  good  for  the  Itmited 
experimental  data  that  are  available. 

We  now  move  on  to  consider  time-dependent  properties  in  the  Hydrogen 
atom-Hydrogen  atom  collision.  Again  we  use  the  same  collision  conditions  for 
the  examination  of  time-dependent  properties  as  in  the  preceding  two  chapters. 
In  this  case,  we  have  a 1 keV  projectile  Hydrogen  atom  colliding  with  a target 
Hydrogen  at  a an  impact  parameter  of  0.2  au.  Figure  10-1  through  Figure  10-8 
show  the  time  evolution  of  the  electronic  populations  in  spin  state  alpha  in  this 
collision  system.  Populations  settle  into  constant  values  for  large  times.  Here  the 
target  atom  is  taken  to  initially  have  an  electron  of  alpha  spin  and  the  projectile 
is  taken  to  have  an  electron  in  spin  state  beta.  During  the  interaction  the  electron 
spins  may  get  exchanged,  and  the  projectile  has  a finite  probability  for  flipping 
its  spin  state  form  beta  to  alpha.  We  notice  the  characteristic  oscillatory  behavior 
of  the  electronic  populations  as  a function  of  time.  In  this  case  the  excitation  is 
predominantly  to  the  2p*  state.  The  beta  spin  electronic  populations  have  very 


similar  structure. 
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CHAPTER  11 

CONCLUSIONS  AND  DISCUSSION 

The  goal  of  this  study  has  been  to  develop  a methodology  to  explore  in 
detail  the  dynamics  of  electron  transfer  in  ion-atom  and  atom-atom  collisions. 
Time-dependent  quantities  describing  these  phenomena  have  been  of  particular 
interest  for  gaining  detailed  insight.  A formalism  has  been  developed  within 
time-dependent  Hartree-Fock  theory  for  the  propagation  of  electronic  density 
matrices  through  time  while  coupling  the  electronic  and  atomic  core  degrees  of 
freedom.  The  eikonal  representation  has  been  used  to  develop  the  form  of  the 
interaction  potential  between  colliding  atoms.  An  averaged  effective  potential 
has  been  introduced  to  ensure  time  reversal  invariance  of  the  potential.  A 
temporal  linearization  procedure  was  shown  to  provide  a method  for  solution  of 
the  electronic  density  matrix  equations  and  to  ensure  the  accuracy  of  this  solution. 
Time-dependent  electronic  populations  and  polarization  parameters  have  been 
introduced  and  calculated  to  allow  for  the  detailed  description  of  the  collisional 
process.  Electron  translation  factors  have  been  included  in  the  calculation  of 
electron  integrals  to  properly  describe  fast  collisions,  and  the  effect  of  the  choice 
of  basis  set  has  been  assessed. 

A new  computer  code  in  Fortran??  has  been  written  to  efficiently  calculate  all 
the  necessary  components  appearing  in  the  integration  of  the  equations  of  the  time- 
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dependent  density  matrix  formalism.  The  flow  diagram  for  the  time-dependent 
density  matrix  computer  code  is  presented  in  Appendix  C.  The  time-dependent 
density  matrix  program  package  has  provided  calculations  on  three  systems  to 
determine  the  validity  of  the  technique  as  well  as  to  investigate  the  behavior 
of  the  time-dependent  properties  we  have  developed.  Calculation  of  charge 
transfer  integral  cross  sections,  state-to-state  integral  cross  sections,  differential 
cross  sections,  and  integral  alignments  in  the  proton-Hydrogen  atom  system  over 
three  orders  of  magnitude  in  projectile  energy  have  shown  a high  degree  of 
agreement  between  the  present  theoretical  results  and  those  of  experiment  and 
other  theories.  The  calculation  of  time-dependent  electronic  populations  and 
polarization  parameters  for  this  system  have  provided  a unique  insight  into  its 
detailed  mechanism  of  electron  transfer. 

Calculation  of  charge  transfer  integral  cross  section  and  state-to-state  inte- 
gral cross  section  in  the  Helium(2-i-)-Hydrogen  atom  system  have  shown  good 
agreement  with  experiment  and  other  theories  for  projectile  energies  at  which 
comparisons  are  available.  Again,  the  calculation  of  time-dependent  electronic 
populations  and  polarization  parameters  have  provided  unique  insights. 

The  Hydrogen  atom-Hydrogen  atom  system,  with  two  active  electron,  has 
also  been  investigated.  The  calculation  of  state-to-state  integral  excitation  cross 
sections  has  shown  good  agreement  with  available  experimental  data.  The 
incorporation  of  two-electron  integrals  in  the  time  propagation  extends  the  scope  of 
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the  application  of  the  time-dependent  density  matrix  code  to  many  active  electrons. 
Time-dependent  properties  in  this  system  again  have  provided  unique  insight  into 
the  collisional  electron  transfer  process. 

The  formalism  of  time-dependent  density  matrix  calculations  is  seen  to  admit 
to  a straightforward  interpretation  of  collision  phenomena.  The  separation  of  fast 
and  slow  degrees  of  freedom  allows  for  a convenient  scheme  for  optimizing  the 
efficiency  of  the  time  propagation  algorithm.  The  ability  to  propagate  the  density 
matrix  through  time  with  varying  time  steps  allows  the  computation  to  concentrate 
more  time  steps  in  the  region  of  the  collision  in  which  most  of  the  change  in  the 
electron  density  is  occurring.  The  definition  of  a reference  density  allows  for  the 
determination  of  the  amount  of  change  occurring  over  any  given  time  step  by 
separately  treating  the  oscillations  in  electronic  densities  due  to  relaxation  of  the 
system  for  fixed  nuclear  positions.  The  explicit  elimination  of  gradient  coupling 
due  to  the  motion  of  the  cores  by  means  of  electron  translation  factors,  assures 
that  the  electronic  population  for  each  state  goes  asymptotically  to  a constant 
value.  This  property  of  our  approach  permits  a detailed  assessment  of  the  role  of 
electron  translation  factors  as  a function  of  projectile  energy. 

Conclusions  can  be  reached  based  on  the  research  presented  here,  pertaining  to 
the  role  of  the  electron  translation  factor  in  the  calculation  of  state-to-state  integral 
cross  sections  and  integral  alignments.  Our  research  indicates  that  ETFs  may  be 
neglected  in  the  calculation  of  state-to-state  integral  cross  sections  provided  the 


110 


core  velocities  are  below  0.02  au  and  provided  the  spurious  asymptotic  couplings 
are  cancelled  analytically.  Further,  our  research  shows  that  the  effect  of  the 
inclusion  of  ETF’s  is  most  strongly  noticed  in  the  integral  alignment  at  higher 
core  velocities.  We  are  able  to  conclude  that  the  differential  cross  section  provides 
the  most  demanding  test  of  the  accuracy  of  the  choice  of  core  trajectory.  It  has 
also  been  found  necessary  to  use  the  more  accurate  averaged  effective  potential 
for  collisions  with  appreciable  integral  cross  sections  at  lower  projectile  energies. 

The  time-dependent  analyses  of  collisional  systems  in  terms  of  electronic 
populations  have  led  to  a number  of  interesting  conclusions.  It  is  noted  that 
electronic  populations  tend  to  fluctuate  during  collisions  and  to  settle  at  a final 
value  only  after  what  can  be  numerous  fluctuations  between  the  colliding  nuclei. 
The  final  value  of  the  electronic  population  is  a function  of  the  projectile  energy 
and  the  impact  parameter.  A change  in  impact  parameter  of  only  0.2  au  can 
make  the  difference  between  complete  transfer  and  no  transfer  in  the  proton- 
Hydrogen  atom  system  for  a 10  eV  projectile.  It  is  seen  that  contrary  to  the 
concept  of  orbital  locking,  introduced  for  neutral  atoms,  the  alignment  of  the 
electronic  system  also  fluctuates  during  ion-atom  collisions.  In  summary,  this 
study  concludes  that  electronic  transitions  in  ion-atom  and  atom-atom  collisions 
are  not  sudden  nor  are  they  necessarily  one-step  processes.  Rather,  electronic 
transitions  exhibit  oscillatory  behavior  and  sometimes  quite  rapid  fluctuations  in 


these  collisions. 
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The  time-dependent  density  matrix  can  be  improved  by  the  inclusion  of  a 
variationally  determined  transition  functional  calculated  in  the  atomic  orbital  basis 
in  which  the  density  matrices  are  written  [30].  Enlarged  basis  sets  including  d- 
orbitals  can  readily  be  incorporated  in  the  time-dependent  density  matrix  code. 
The  description  of  nuclear  motions  can  be  generalized  in  practice  to  include  core 
electrons  as  well  as  nuclei.  Finally,  the  method  is  suitable  for  the  inclusion 
of  more  than  two  cores  and  could  be  used  for  collision  phenomena  involving 
diatomics  and  polyatomics. 


APPENDIX  A 

CALCULATION  OF  OVERLAP  AND  HAMILTONIAN 
MATRDC  ELEMENTS  IN  THE  HYDROGENIC  BASIS 

In  this  appendix  we  write  the  two-center  overlap  and  one-  and  two-center 
Hamiltonian  matrix  elements  for  a Hydrogenic  basis  including  Is,  2s,  2pi,  2po, 
and  2p_i  orbitals  on  each  center.  These  spherical  representations  of  the  p-type 
functions  are  related  to  the  Cartesian  in  the  following  way 


Integrations  are  carried  out  in  elliptical  coordinates  as  defined  in  chapter  7.  The 
centers  are  referred  to  as  center  a and  center  b and  orbitals  are  referred  to  by 
labels  with  subscripts  to  denote  the  center  on  which  the  orbital  resides.  Hence 


represents  the  overlap  between  any  orbital  on  center  a and  the  Is  orbital  of  center 
b.  The  expressions  for  integrals  are  tabulated  as  functions  of  the  intemuclear 
separation  R and  the  effective  charge  of  the  centers  Cm-  The  following  integrals 
are  introduced  in  the  interest  of  compact  notation. 


I2pi)  = 

|2py)  = - 12p-l)) 

12pz)  = l2po) 


(A.1) 


(X'oll'St) 


(A.2) 


(A.3) 
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I 

B(/3)  = J T]exp{-0ri)dT] 


-1 


= (-^  - ^)exp(-^)  + (-i  + ^)exp(^) 


Id 


i 

C{0)  = J exp{-pT])dT] 


-1 


i 

D{ld)  = j T]^exp{-0ri)dT] 


-1 


. 1 3 6 6 ^ 

, 1 3 6 6,  ... 

1 

■P’(/^)  = J »?^exp(-/?r/)d7/ 

4 12  24  24 


-1 


= (■ 


> /?2  13^  13^  15^ 


)exp(-^) 


/I  4 12  24  24.  ... 

G{a)  = J exp{-a()di  = ^ 


OO 

//(a)  = J iexp{-a^)d^ 


= (-  + ^)exp(-Q) 

Q Q' 


(A.4) 


(A.5) 


(A.6) 


(A.7) 


(A.8) 


(A.9) 


114 


00 


J{a)  = j i^exp{-a(]di 
1 

/I  2,2-  . . 

= (-  + — + —)exp(-Q) 


a Of 


00 

K{a)  = J e\p{-a^)di 

1 

A 3 6,6. 

= (-  + -2  + -3  + — 
Q a-* 

00 

L{q)  = j exp(-aO<^^ 


(A.IO) 


(A.ll) 


(A.12) 


1 4 12  24  , 24.  . - 

= - + ^ + ^ + -4  + — )exp(-Q) 
a O'*  Q® 

The  following  is  the  list  of  overlap  integrals  in  terms  of  the  above  integral.  Note 


that  in  general 


(Xa|\6)  = (X6lXa)  (A.13) 

as  the  hydrogenic  orbitals  are  real,  so  only  the  first  of  the  pair  of  forms  will  be 
reported.  Also  if  the  two  orbitals  are  on  the  same  center  they  are  orthonormal. 

{ls.|ls»)  = [/!(/?)  J(a)  - CWG(a)] 

(l»,|2s,)  = - C{l3)G{o)] 

- - CWH{a)  + W)G(o)| 

o = j(2c. + a)R  i e = j(2<.  - m 


(A.15) 
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(U.|2to)  = - C{H)a{a)  - B{l3)I<W  + DWH(o)] 

a = ~i^Ca  + Cb)f^  » ^ — T(^Ca  ~ Cb)^ 

4 4 (^16) 


(23,|2si)  = (Wy(a)  - C(^)G(a)]  - ^[A{ff)KM 

+ B{^)J{a)  - C{I3)H{q)  - D{^)G{a)]  - ^[A{0)K{a)  - B{l3)J{a) 

- C{P)H{q)  + D{l3)G{a)]  + - 2C(/3)J(a)  + F(/?)G(a)]} 

Q = T(Ca  + Cb)B  ; 0 = jlCa  ~ Cb)B 
4 4 


(2».|2pot)  = (l^(^)J(o)  - C(/3)G(a)  - B(3)A'(q) 

+ D{P)H{a)]  - ^[/l(/S)A(o)  + B(/3)J(a)  - C()J)A(a) 

- Z)(/J)G(q)  - B(/3)i(a)  - C(^)A(q)  + D(y3)J(a) 

+ W)A(a)l) 

o = 7(C.+a)Ai  ^=7(C.-C»)A 

4 4 

{2pia\2pib)  = - C(^)I(a)  - A(^)J(a) 

+ F{l3)J{a)  + C{0)G{a)  - F(y3)G(a)] 

a = i(Ca  + a)^;  /3  = ^(C«-a)^ 


(A.18) 


(A.19) 
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(2po.|2po»)  = - CWGM  - CmM  + F{P)J{a)] 

o = t(C.  + <b)F ; ^ = t(C.  - (b)R 

4 4 (^20) 

All  other  overlap  integrals  are  zero  from  integration  over  the  y?  variable.  The 
following  list  contains  the  non-zero  one-electron  Hamiltonian  matrix  elements. 

{Ua\H\Ub)  = Ei,,{Ua\lst)  - - B{0)G{a)]  ^i) 

a = ^(Ca  + Cfc)^;  ^ = \{Ca-Cb)R 


\/C°C^  R~ 

{2s^\H\Ut)  = {2».|lsj)£',„  - ^ {AWH(,a) 

- B{/3)G(a)  - ^-i^\A{0)J{a)  - C(/J)G(a)]} 

a = j(C.  + i<b)R ; e = j(C<.  - Kb)R 


(Ml) 


{2poa\B\^^b)  — (2p0a  |l^fc)^lst  16^ 

- B{l3)G{a)  B{^)J{o)  - C{i3)H{a)] 

a = l(c«  + i<b)R ; /9  = i(C.  - ib)R 

4 4 


W)H{a) 


(A.23) 
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{2Sa\H\2Sb)  = {2Sa\2Sb)E2s,  - 


{A{l3)H{a) 


- B0)G(a)  - ^[A(B)J(a)  - C(0)G(a)]  - ^lA(ff)J(a) 

- 2B(P)H(a)  + C(^)G(a)l  + ^i^(A(^)A-(o)  - B(B)J(o) 

- C(«A(a)  + D(/9)G(a)l) 


a 


= j(Ca + m \ 13 = j(Ca  - m 

{2poa\H\2sb)  = {2poa\2sh)E2sb  - 


{A{l3)H{a) 


- B{B)G(a)  + BW)J(o)  - CWH{a)  - ^[A(S)J{a) 

- 2B{i3)H{a)  + C(/3)G(q)  + B(/3)A(a)  - 2C(,9)J(q) 


+ D{B)H{o)]] 


a = + m ; » = j(C.  - &)Ji 

(2poo!^|2pot)  = (2p0a|2po6)£'2po6  “ 


W)H{a) 


- 5(/3)G(q)  - C(^)A'(q)  + D[^)J{o)] 

o = i(c„  + a)^;  l^  = kCa-Cb)R 

4 4 

ICC^R^ 

{2pu\H\2pa)  = {2pu\2pu)E2p„  - |B(/J)G(a) 

- A{0)H{a)  + A(B)K(a)  - B(B)J(a)  + C(^)//(a)  - D(B)G{a) 

- Cm<{a)  + C(;9)J(o)] 

o = 7(C.  + (t)R ; e = 7«.  - Ci)fl 

4 4 


(A.24) 


(A.25) 


(A.26) 


(A.27) 
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{U.\H\U.)  = Eu.  - S<^[A(i3)H(a)  + B(^)C(a)] 
a = 0 = Co-R 

(2s,|fl|lsi)  = + B0)G{a) 

- ^W0)JM  + W(0)H(a)  + C{ff)G{a)]] 

a = ^ = ^CaR 
4 

(2po.l«|l^a)  = -^^I/l(/3)//(a)  + B(^)G(a) 
BWJ(o)  + C(/3)K(a)l 

c.  = 0 = \uR 

4 

(25a|^l2s„)  = - ^^{A(^)//(a)  + D{B)G{a) 

- ^[AmJ{cx)  + + C(^)G’(a)] 

+ (a)  + W{^)J{a)  + 3C(^)//(q)  + D{0)G{a)]] 

16 

a = ^ = ica/? 

(2.a|^|2poa)  = -^^{A{0)H{a)  + B{B)G{a) 

+ B{l3)J{a)  + C{B)H{a)  - ^[A{B)J{a)  + 2B{l3)H{a) 

+ C{B)G{a)  + B{l3)K{a)  + 2C{0)J{a)  + D{l3)H{a)]} 
a = B=  \CaR 


(A.28) 


(A.29) 


(A.30) 


(A.31) 


(A.32) 
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(2p„.lW|2po.)  = E2,,  - ^^^\A(l3)H(a)+  B(»)G(a) 
2B(P)J(a)  + 2C0)H(a)  + C(/3)A'{q)  + D{0)J{a)\ 
a = ff  = ^(,R 

(2p,.|A|2p„)  = E.,.  - S^[-A{/3)H(o)  - BWCM 
+ A(ff)K{a)  + B{/3)7(q)  + C(e)H(a)  + D(»)G{a) 

- CWK(a)  - DmJ{a)] 
a = 13  = jCoB 


(A.33) 


(A.34) 


APPENDIX  B 

CALCULATION  OF  ONE-ELECTRON  INTEGRALS  OVER 
GAUSSIANS  INCLUDING  ELECTRON  TRANSLATION  FACTORS 

In  this  appendix  we  present  the  two-center  overlap  matrix  elements  and  the 
two-center  one-electron  Hamiltonian  matrix  elements  in  a basis  of  Gaussian  or- 
bitals including  electron  translation  factors.  The  general  problem  under  consider- 
ation is  the  evaluation  of  an  integral  of  the  form 


We  consider  the  case  in  which  f and  g are  real  Gaussian  functions  and  exp{iq  • r) 
corresponds  to  an  ETF.  The  integration  is  earned  out  in  Cartesian  coordinates  and 
Gaussians  of  differing  angular  momenta  are  related  to  each  other  by  differentiation 
with  respect  to  the  coordinates  of  their  centers.  We  consider  s-,  p-,  and  d-type 
Gaussians  which  are  defined  by  the  following  relations 


(B.l) 


(B.2) 
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As  discussed  in  chapter  7,  we  need  the  two-center  overlap  integrals  which 

h,.  = )exp{iq-  rj)gs{l3,f„j) 

= (-)t(a/?)</i  exp(— j-  + iq-  Rp) 

K 

. n n ^Rm  "I"  ^ Rn 

K = Oi  p Rp  — ') 

K 

K = exp[(— - -Rn)‘] 

K 

= j ^j9pi  ( a , Tm J ) exp ( uy  • Vj )gs{^,  rnj ) 


'P<.» 


,2q^ 


— )(Cm  Cn)  + ]^s,s 

\yOl  K h. 

h,pt  = J (t‘rjg,{a,fm,)exp{iq-f,)g,^(l),fpj) 

= -L|(M)(^„  _ c.)  + iM]/, , 

VP  K K 

^PoPr,  = J d^rjgp^{a,fmj)exp{iq-fj)gp^{l3,f„j) 

1 2o^ 

= ~1=^\ — 

Vap  K 

+ [ ((.m  ~ ^»t)  "I  ][  (7m  7ii)"h  , Jj-^s.s 

/C  K K K 


hir,,s  = J d^ngd^,^{a,frn])exp{iq-fj)gs{l3,r,,j 


) 


- i/tzMu 

+ [-(  — )(Cm  - Cn)  + ^][-(  — )(7m  ~ 7n)  + 

K K K K 


Js,d(r,  = J d^rjgs{a,frrij)expiiq-fj)gd^^{l3,i^„j 


-inZMu 

" 13^^  K 

+ [( )(Cm  - Cn)  H ][(~“)(7m  “ Vn)  H — ]}/5,a 

/C  /C  /v  ^ 


are 


(B.3) 


(B.4) 


(B.5) 


(B.6) 


(B.7) 


(B.8) 
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= j d\j9d^^{a,f^j)exp{iq-rj)gp^{0,fnj) 


+ 


2ia^ 


(aq^S^p  ^qt]^(p)] 


,2q^ 


+ {(-^)<5c.  + [ 


— 4q^/9^ 


ArA 


C^v 


ajdqcqrj  , 2ia^ 


+ ■ 


[^A,  + i2££|)/.,, 


A(r  — (^m  ~ ^n 


^P(Ap  = j d^^j9pi  (a,  r„,j ) exp(rg  • Tj )9d^^{P,  r^j ) 
{[— f(6cpA,  + Ac6,p) 


+ 


^y/a 

^^{-aq^Srip  + /?9,<5cp)] 


,2Q;9 


+ {(~~)<5cv  + [ 


«2 


Ac  At, 


a^qc,qr)  ‘2ial3 


+ 


K? 


■(q9(^A, 


K 


Ap  + 


/^9»?Ac)]}  X 


(B.9) 


^9»jAJ]}x 


(B.IO) 
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h,„dr,a  = J d^rjgd^^{a,fmj)exp{iq-fj)gd,Al^,f„j) 

1 3'^ 

= ^{^(<5cpV  + M.p) 

ap 


.4a^^2  2ia3.  ^ „ , „-2a/3  .iqp^ 


2a^ , 

—kr,  - 

-Aa^3^  A A . 2m^ 


2a/3  4a2^2  q^9c9^ 

+ i ^ ®C»/  ^2  ^2 

+ ^(a,(A,  - 5„Ac)][^A,  + i^])  X 

^ /C 


(B.ll) 


As  discussed  in  chapter  7,  we  also  need  the  two-center  one-electron  kinetic  energy 
operator  matrix  elements  in  the  Gaussian  basis.  These  elements  are  as  follows 


Ts,s  = [3a  - 


3a*  2a^3‘ 


— (ii™  - iin?  + U 


9 


(B.12) 


3a  2a2/?2  - 


^p<.«  ~ 


2^2 


{Rm  - Rn?  + ^ 


2ia~3 


+ — ^ — q'  {Rm  ~ Rn)]Jpi,s 

3 

2ia2  3q^ 


+ [ 


—iai  3^ 


K- 


(Cm  Cn)  T 


o 


]/s,a 


(B.13) 


^s,p<  — — 

K 

2m^/3 


2fc2 


H 2 ^ ' (-^”1  ~ •^«)]'^«,Pc 


+ [ 


Aa^j3'^ 


(Cm  Cn) 


2ia‘^y/^q( 


]l 


s.s 


T ^.^(R  -R)^^ 

^P<,Pr,  — ^ ^2  2k” 

2io^ j3  -,  , ft  f5  \i  r 

H ^2  ^ ' (-^m  ~ •^n)]-^p^,p,, 

/c** 

K*  K 


+ [ 


2ai 


{llm  Vn) 


K? 

3 ^3 
4q2^2 

[ ~2 


2m"\/^9T,i, 

■J-'Pc>« 


0 

K- 


Td.„,s  = [3a  - 


3a2  2q2^^ 


2„2 


+ 2^,--  (4  - 4))/j„,, 


2^2 


, r ■4a^^^„  ..  \ , 

+ [ 2 ..2  J‘'P<:>^ 


(B.14) 


(B.15) 
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3a“  ^ 


0 0 

2 . 


2/c2 


+ [ 
+ [ 


2/3# 


4q^/9" 


•(j/m  ^n) 

"(Cm  ~ Cn)  “ 


«2 


U.,: 


Pc 


2ia-y/^q(; 


0 VSTH  STl/  9 

/C^  K" 


Pv 


-Aa-I3  , 

+ [ — “2 — kv]^s,s 


T - rq  _ ^ _ ^ _ o , i_2_ 

^d<P,Pr,  ~ ^ ^2  '*’  2«2 


2 . 


+ 


2ia^^  _ f! 


«2 


q-  {Rm  Rn)]Id^p,p,, 


+ [- 


-4q2 


(/^m  /^n)  “t* 


2m2/?<jf^ 


9 vrm  r7t;  ' 9 

K- 


l^K,: 


Pc  .P'7 


+ I^(C.,-C,.)  + ^1/.., 


K- 

Aa-0l 


2ior^f^qr^ 


+ r-(^/m  - J?n)  - ,.2"  *^]-^^CP.^ 

K f\ 

-Aa0^  , 

+ I — ~2 — ^Cp]'^SP'7 

fX 

4Q2/52 

■*■  I ^.2  °whp.s 

fi 

^4q2/9§  ^ ^ ^ 

+ 1 ;32  ^’>pj-'pc>« 


(B.17) 
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~ ^ 2«2 
- ^q-  (Rm  - R,)Vn4,. 

, ,4a2/?t^  2iory/^qp 

+ 2 — ~ ^"'  1^2  i-'PC.P'i 

/C^  ^ 


+ 1 


4q^/32 


(^m  ^n) 


2ia~\/^qjj 


9 V7TTI  7W/  9 

K“^  K“ 

Mt 


Kk.- 


Pi  tpp 


Aa^/3'^ 

+ [ — “2 — ^Cp]-^«p-7 
t\ 

3 3 

Aq^P^  , 

+ I “2  ®C»/J^«.Pp 

lx 

+ [ “2  ®»?pJ-'Pc-s 


(B.19) 
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Tdi^dr,.  = ~ 2ac2 


+ 


2ia^j3 


Q • (^m  -Rn)]-^d,-p,(f 


^pi“T7<f 


-4a2/3^ 


r , X 2mi^9p, 

+ [ ;;;2  «.2  J-'Pc.'^-i 


+ [^(U-C»)  + ^1Va. 

, r4a2/9§,  , 2ia-y/^qr,^j 

I ^2  ,c2  y<^(P<P<r 


+ [ 


2/34 


■(<^m  ^n) 


2io“s/^qa 


ipiPv 


-4a/3“ 

+ [ — ~2 — kpVs.d„a 

fV 

3 3 

^4a2/^2  ^ ^ ^ 

+ [ 2 °p<^\hi,Pf} 

rx 

3 3 

.4Q2^2 

+ I ®C'?J-'Pp,P<t 

lx 

3 .3 

^402^2  ^ ^ J 

+ I ^2  ®ppi''pc.p<’ 

^4a2/?2  ^ ^ ^ 

i ^2  ®C<^J-'Pp.Pti 

tx 

-Aa-/3  , 

+ 1 ;;:2 


(B.20) 


APPENDIX  C 

FLOW  CHART  OF  THE  TIME-DEPENDENT 
DENSITY  MATRDC  COMPUTER  PROGRAM 


Advance  time  to  t=tO+dt=tl 
Find  positions,  velocities,  and 
PO(tl) 


Calculate  H,  S in  body-fixed 
frame  and  rotate  to  space- 
fixed  frame 


No 
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